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HODELING CHARGE DIFFUSION FROM PROTON- INDUCED IONIZATION 
TRACKS IN SILICON PHOTODIODE ARRA/S 

By 

L. Cole Smith^ and Jacob Bechet^ 

ABSTRACT 

We are modeling Che combined problem of generation and collection of 
charge carriers created during passage of energetic protons through a sili- 
con photodiode array. Vife have also experimentally obtained pulse-height 
distributions of noise charge collected during exposure of a digicon-type 
diode array Co 21 and 75 MeV protons, It is shoyn Chat the magnitude of 
charge collected by a diode from each proton event is determined not only by 
diffusion, but by statistical considerations involving Che ionization 
process itself. Utilizing analytical solutions to the diffusion equation 
for transport of minority carriers, together with the Vavilov theory of 
energy- loss fluctuations in thin absorbers, we present simulations of the 
pulse-height spectra which follow the experimental distributions fairly 
well, and provide us with an estimate for the minority carrier diffusion 
length Lj. 


INTRODUCTION 

The introduction of LSI microcircuit technology has made possible Che 
construction of planar arrays of small p-n junctions on the surface of a 
single thin chip of silicon. Since the depletion zone of each junction can 
behave either as a tiny capacitor, or as a potential well, these arrays can 
be used to store patterns of charge which can be retrieved quickly and, if 
necessary, repeatedly. These devices are conveniently divided among two 
categories: (1) analogue, including image sensors (such as linear photodi- 

ode arrays (LPDA) and charge coupled devices (CCD's), and (2) digital, the 

^Graduate Research Assistant, Department of Physics, Old Dominion Univer- 
sity, Norfolk, Virginia 23508. 

^Associate Professor, Department of Physics, Old Dominion University, 
Norfolk, Virginia 23508. 


prime example being semiconductor memory chips (e^g. MOS dynamic RAM's (rair* 
dom access memory). Advances in LSI and VLSI techniques have increased the 
number of diode elements that can be arranged on a chip, by decreasing the 
size of each element, so that 64- kbit dynamic RAM's are now found with 
individual memory cells of only a few square microns in area. Although 
smaller device dimensions are likely to be forthcoming as the technology 
improves, it has been demonstrat&e in recent years that as downscaling 
continues, most of these microelectronic devices b?>corae critically sensitive 
to energetic particle radiation, such as found in cosmic ray showers, and 
natural radioactaivity on earth and the Van Allen radiation belts and 
primary cosmic ray particles above the atmosphere. This consideration 
represents a very serious operational problem associated with the use of 
state-of-the-art VLSI technology in both analogue (imaging) and digital 
(memory) circuits, for either ground or space applications. 

Radiation Effects on Digital Devices 

Considerable attention has been paid recently to the problem of soft 
error generation in high-density dynamic RAM chips (refs. 1-5). (By soft 
error we mean random, nonrecurring, single bit errors in memory devices. 

The errors are not permanent). May and Woods (ref. 1) and Yaney, Nelson and 
Vanskike (ref. 3) demonstrated that alpha particles emitted during the spon- 
taneous fission of U 238 and Th 230, (present in parts-per million levels in 
the ceramic packaging materials typically used for IC's), are capable of 
generating enough electron-hole pairs in the depletion region and silicon 
bulk to discharge a single cell on an LSI chip. Their results showed, pre- 
dictably, that the alpha-particle- induced error rate depends upon the cell 
size, the energy of the alpha particle, the charge collection efficiency of 

each cell, and the critical charge Q i.e., the number of stored elee- 

xrit 

trons which differentiates between a "1" and "0" (ref. 3). The critical 
charge for the 16 K RAM’s used in the experiments was estimated to be about 
200 fC, which required an alpha particle having an energy of at least 4.5 
MeV to cause an upset. Other experiments have shown that CCD memory arrays 
are even more sensitive to ionizing radiation because, due to the increased 
levels of integration accompanied by a reduction in storage cell size, they 
possess a smaller value for (ref. 6). Beta and gamma sources produce 
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negligible loft errors in dynamic memories because of their relatively lower 
ionization loss rate. Energetic ions may interact through a (p, alpha) or 
similar nuclear reaction to produce soft errors, but this effect decreases 
rapidly with decreasing particle kinetic energy (refs. 7, 8). Cosmic-raj^ 
induced latch-up in CMOS logic and memory modul/j!B appear to remain a serious 
operational problem since its discovery onboard Air Force satellites in 1975 
(ref. 9) . 


Radiation Effects on Analog Devices 

The effects of particle radiation on image sensors has received a less 
exhaustive treatment, possibly due to the more widspread use of the ubiqui- 
tous memory chip. Multidiode image arrays are still fairly expensive, and 
users of these devices are perhapp unwilling to subject them to radiation 
levels that could be potentially destructive. Actually, tests have shown 
image sensors to be equally sensitive to extremely low fluences of particle 
radiation (refs. 6, 10, 11) . Also, unlike a memory cell, which experiences 
only a change of state, each diode element is, by design, sensitive to a 
wide range of charge levels, down to about 0.5 fC. As a result, its sensi- 
tivity to charge noise is not limited to the effects of alpha particles, but 
extends to nearly any charged particle losing sufficient kinetic energy to 
produce ionization in the silicon. 

Modela of Charge Collection in Integrated Circuits 

All of these considerations have inspired several authors to propose 
various charge diffusion models to explain the various effects resulting 
from ionizing radiation passing through silicon microelectronic devices 
(refs. 12-15). Techniques range from simple analytical solutions to the 
diffusion equations (12-15) to full-scale Monte Carlo simulations (refs. 13, 
14), and so far results appear encouraging for predicting the soft error 
rate due to alpha- particle penetration of the chip. 

To our knowledge no previous theoretical work has treated the problem 
of charge collection from proton-created ionization tracks in integrated 
circuits designed for analogue charge detection. Part of the problem is 
that generally little is known experimentally about charge collection from 
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particle tracUa in these structures* Consequently several approximationa 
comnonly used in alpha particle work have gone untested* For examplei many 
authors assume minority carrier recombination plays little part in dk^termin* 
ing the amount of charge reaching a diode element from the path of a charged 
particle. Kirkpatrick (ref. 12) as 'unes recombination lifetimes of several 
milliseconds, although experiments with 4.2 MeV alpha particles bombarding 
CCD arrays have indicated minority carrier lifetimes t of leas than a 
microsecond (ref* 6)* This corresponds to a diffusion length of 26 U, using 
Lj "yDT, which still exceeds the range of most alpha particles in silicon. 
Protons# however, can produce tracks equaling the thickness of the silicon 
wafer (100-400 m) since they essentially penetrate the chip with little 
energy loss. This possibility of charge diffusing to the surface from very 
large distances necessitates the inclusio;i of recombination effects in 
charge collection models involving protons, 

Another consideration related to proton-induced ionization in silicon 
involves energy-loss fluctuations along the ionization track* This effect 
is significant only for particles which deposit a small fraction of their 
kinetic energy in the medium, and hence has not been treated In alpha-parti- 
cle models, However, energy-loss straggling plays a major role in determin- 

ing the maximum amount of charge collected by a diode from proton-induced 
ionization tracks* 

Experimental Work with the Digicon Diode Array 

Both CCD arrays and LPD arrays are used as photoelectron targets in 
many image tube camera systems developed over the past few years. Some of 
these cameras, such as the Digicon, are to be flown on space astronomy 
satellites in the mddle 1980' s, This will expose the imaging components 
(including the diode arrays) to energetic particle radiation potentially 
capable of causing an interfering noise background. So far no arrays have 
flown in space, so little in situ data are available. Therefore, this 
information has to be obtalited on the ground, using man-made particle beams 

of appropriate energy to simulate the radiation environment. 

We performed these tests for NASA/ Goddard in 1979-80, vdien three 
different 512 element linear photodiode array were irradiated with proton 
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beams at the University of Maryland 100 MeV If^'^hronous Cyclotron, Two of 
those arrays were the targets of engineering model Oigicon tubes, similar to 
flight models to be flown on the Space Telescope, The other array was part 
of a breadboard circuit but was similar to the other two, The purpose of 
the experimental work was to evaluate the effect of South Atlantic Anomaly 
levels of proton radiation on the performance of the tube. This paper, 
which includes data from exposures to protons at 11,5 and 75,4 MeV, uses the 
essential results of );ho8e tests as a basis fur the evaluation of a diffu- 
sion model that includes the effect of the shape and pattern of the diodes 
on the array, mirority carrier recombination, and finite recombination 
velocities at the backside of the silicon wafer. Three-dimensional simula- 
tio'ns are used to shown the necessity of including the effects of energy- 
loss straggling in accounting for the experimetal results. The intent of 
this model is twofold: first, to provide a solid theoretical framework 

within which to evaluate p«svticle radiation effects in multi-diode struc- 
tures, and second, to use this framework to investigate the dependence of 
experimentally measurable quantities on physical parameters in silicon. 

Wo have emphasized the negative aspects of the charge sensitivity of 
analog silicon devices. It should be pointed out that the charge sensitive 
nature of the diode array coupled with the precise geometrical arrangement 
of diode elements could prove useful in nuclear and solid state experiments, 
with the possibilities of position sensitive detection of high energy parti- 
cles and studies of carrier transport phenomena in silicon. Considerable 
attention has been devoted recenti , to the study of the determination of 
transport parameters of minority carriers in single p-n junctions by employ- 
ing scanning electron beam excitation. Calculations have appeared in recent 
years in association with electron beam induced current (EBIC) experiments 
designed for quantitative determination of diffusion lengths, bulk lifetime 
and surface recombination velocity of minority carriers in silicon (refs, 

28, 29, 30). However, none of these results have been applied to under- 
standing charge diffusion and collection in microelectronic chips involving 
many diode elements. The extension of this work to multi-diode arrays ap- 
pears promising. As recent tests utilizing cosmic-ray showers have indi- 
cated, the CCD array, because of its charge imaging property, appears to be 
especially suitable for experiments measuring carrier diffusion lengths and 
surface recombination effects (ref. 7). 


Hence, both a theovetlcal and experi>nental void needi to be filled In 
order to understand the effects of Ionizing radiation on these devices* 




THE SILICON DIODE ARRAY TARGET 
Introduction 

Before we proceed with «n easei^nent of the reiponee of photodiode 
erreyi to particle redietioni a brief deacription of the fabrication and 
atructure of theae devicea ia given. Particular enphaaia ia placed on 
arraya aiuilar to thoae uaed in our experinMintal work, 

Ik 

Diode Array Fabrication 

The enploynent of ailicon-dioder*array targcta in varioua electron- 
imaging modea of operation, particularly lot^*light level applicationa, haa 
been diacuased in the publiahed literature (refa. 16, 17, 18), Moat multi- 
diode atructurea utilize a one- or two-dimenaional array of aenaing elementa 
in which (for the caae of a aolid atate image aenaor) incident photona or 
photoelectrona are abaorbed and converted to electric charge. The aenaora 
or atorage elementa may be diffuaed photodiodea or MOS capacitors (surface 
depletion layers induced by a voltage applied to semi-transparent electrodes 
of a material such as polycrystalline silicon). The latter corresponds to 
the conventional CCD sensor array. These sensing eleinents are created on 
substrates made from thin slices of either n or p type silicon, usually 
homogeneous, (although epitaxial wafers are nowadays not uncommon). Bulk 
resistivity is usually low, around 10 ohm-cm. Since we are considering only 
photodiode arraya, we shall limit our discussion to their fabrication. 

In preparation for forming diodes, the wafer surface is thermally oxi- 
dized ('^1,000* C) . Diode elements are formed using standard photolitho- 
graphic and planar processing techniques employed in the manufacture of 
integrated circuits. This includes etching away the oxide layer according 
to the desired array pattern and forming p + n junctions by the diffusion of 
boron into the n-type silicon through the open windows. The shallowness of 
these diffusions (typically 1 y), coupled with the low resistivity of the 
substrate creates very thin depletion layers when reverse bias is applied 
to the substrate. For photoelectron detection, only a few microns of 
sensitive volume is necessary so these thin junctions (typically 3-4 U) are 
adequate. The remainder of the substrate serves as a rigid support for the 
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array and as a surf/ice on which to bond the diode signal Leads. Substrate 
thickness can vary from 10 M to nearly 500 M, depending on the application. 

A cross-section of a typical diode array is shown in Figure 1. 

Detection of Charge Signal 

Methods for detecting the charge collected at each diode are almost as 
numerous as the different types of sensors. Analogue and digital shift 
register scanning of the diode elements are used in both CCD and photodiode 
arrays; current generation LPD arrays use serial multiplexing whereby each 
diode is sensed sequentially. Charge Injection and charge transfer tech- 
niques are usually associated with devices using MOS capacitors as sensing 
elements, and will not be discussed here. In photodiode arrays, conversion 
of the charge into a measurable signal in the external circuit is accom- 
plished in one of two ways: charge- storage and charge- integration. 

Charge- storage operation is based on the fact that if a p-n junction is 
reverse biased, then open circuited while illuninated with a light source, 
the charge initially stored on the depletion- layer capacitance decays at a 
rate proportional to the incident illumination level. Thus, by monitoring 
the charge required periodically to re-establish the initial-voltage condi- 
tion, one will obtain a single proportional to the incident illumination. 
This method is not suitable for extremely low levels of illumination, when 
only single photons are incident, as the photon- induced diffusion currents 
are so small as to be insufficient to discharge the charge stored in each 
diode's potential v^ell. In this case some intermediate gain stage is neces- 
sary, usually in the form of a photocathode and an accelerating potential of 
around 20,000 V. 

In the latter case, a photon- ejected electron replaces th*i photon as 
the means of creating diffusion currents in the silicon, an.<i the charge 
integration method is used to detect the signal. In thi? method charge 
sensitive preamplifiers connected to each diode are used to integrate the 
brief current pulse resulting from a single impinging photoelectron. The 
resulting voltage signal is thus proportional to the total charge either 
released in or eventually reaching the diode depletion layer. Using conven- 
tional pulse shaping and counting techniques, the incident illumination 
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level et the photocathode can be determined asauning knowledge of auch 
factora aa quantum efficiency, diode area, etc. 

Radiation Senaitivlty 

Although uaually deaigned to abaorb photoelectrona having energiea 
around 20-50 kiloelectron volta (keV), (depending on the accelerating volt- 
age at the photocathode), each diode in the array will obviously detect any; 
charged particle losing enough energy in the silicon to produce a signal 
exceeding the electronic noise level (typically 4 keV equivalent noise 
charge). In fact, each diode can collect charge from two regions*. (1) 
inside the diode's depletion layer, where large electric fields (10^ V/m) 
cause rapid collection (10 nsec) of charge cteatad there, and (2) outside 
the depletion layer, where there are no electric fields, and charge collec- 
tion occurs through the relatively slow process of minority carrier diffu- 
sion in the homogeneous n (or p) substrate (see Figure 1). As a result, not 
only is an individual diode cell's noise sensitivity enhanced, but on an 
array containing many diode elements spaced closely together, the deeper an 
ionizing particle penetrates into the silicon wafer, the greater the poten- 
tial for activating many diodes simultaneously, That this indeed takes 
place has received little attention in the literature. The author is aware 
of only two papers that experimentally document the effects of particle 
radiation on image sensors. Ko (ref. 6) and Marcus (ref. 11) briefly 
describe the effects of alpha particles and secondary cosmic ray partaicles 
on CCD sensor arrays. Marcus' work with the Fairchild CCD-211 shows vivid 
detail of charge diffusing over dozens of pixels resulting from the ioniza- 
tion tracks of secondary charged particles as they penetrate the wafer 
during the course of a 75-min exposure. Marcus presents no quantitative 
analysis, but Ko uses the results of similar experiments with alpha parti- 
cles to estimate the minority carrier diffusion length of the substrate of 
his CCD array. It appears at present that our work, described next in 
"Experimental Work," represents the only data available on the effects of 
monoenergetic protons on imaging photodiode arrays. 
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EXPERIMENTAL WORK 


In Produce ion 

The silicon diode has been used for many years by nuclear physicists 
and others to detect electrons, protons and heavier particles from radioac- 
tive sources and accelerators. Fast, lof^*noise pre- anplif iers have been 
developed to pei^nit energy determination with high resolution. However, 
surprisingly few data are available on the response to particle radiation of 
semiconductor devices inteujed for other purposes e.g, silicon photodiode 
arrays, CCD's, etc. This despite the ready availability of many of these 
devices ap well as their potential usefulness as an experimental tool for 
the solid state physicist. 

The response of a 512 element linear photodiode array to mediun energy 
proton radiation has been described elsewhere by the author (refs. 19, 20, 
21). The array was the photoelectron target of a Digicon tube similar to 
one being flown on the Space Telescope in 1985. In our test setup, the 
array was irradiated ^ situ by a focused monoenergetic proton beam while 
selected diodes on the array were monitored. Signals induced by the protons 
penetrating the diode array were measured, digitized and counted using pulse 
height and coincidence techniques described in (ref. 19), (see test set-up 
in Figure 3) . 


Description of Digicon Camera 

Although the Digicon has been described in review articles elsewhere 
(refs. 22, 23, 2A) , a brief description here is in order. Essentially a 
digital imaging device, the Digicon is a multichannel parallel- output photon 
detector which uses a thin silicon diode array target to detect accelerated 
photoelectrons from a UV or visible sensitive photocathode. The photocath- 
ode is deposited on the inside surface of a transparent window, in this 
instance, magnesium flouride ()%F 2 ) . In normal operation a spectrum image 
is focused onto the transmission photocathode at the front end of the detec- 
tor, (see Figure 2). Ejected photoelectrons are accelerated through a 20-22 
kV potential and are magnetically focussed by a rare earth permanent magnet 
assembly onto a linear target array of 512 discrete photodiodes. (In our 
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test set-up, this magnet assembly was replaced by a solenoid). Each of 

these diodes consists of a thin p + n junction imbedded on the surface of a 

400 p thick, n-Si (p ■ 4.5 ohm- cm) homogeneous substrate to which a +5 V 

reverse biasing potential is applied. Each diode measures 40 P wide by 250 

p long and is spaced at 50 p intervals. The geometrical arrangement of the 

diodes on the array, as well as the shape of each diode, is shown in Figure 

4. When biased, each diode forms a depletion layer having a depth of about 

2.5 M . Every photoelectron that is stopped in the depletion zone results in 

the creation of a charge pulse of several thousand electron hole pairs. The 

presence of large electric fields in this region results in rapid separation 

and collection of this charge onto the aluminum electrodes deposited over 

about a third of the surface of each diode. Each of the 512 diodes is con- 
» 

nected to a separate channel containing a charge sensitive preamp, shaping 
and amplification circuits and signal- processing electronics which digitize 
the measured charge from each diode independently into a digital word for 
storage in memory. 

Test Configuration 

The. test configuration is shown in Figure 3. Collimated proton beams 
of 1 cii? cross-section were directed onto the center of the faceplate at the 
front of the Digicon, which was located inside a scattering chamber and 
securely bolted onto a movable arm. The Digicon was free to rotate around 
an axis centered on the photocathode, allowing the rear of the tube to be 
swung out of the path of the proton beam. This allowed monitoring of the 
luminescence in the MgF 2 faceplate without the strong interference created 
by protons penetrating the window and striking the diode array. To monitor 
the direct proton bombardment of the array, the arm could be returned to 
zero degrees very quickly, and the high voltage removed. With the photo- 
electron background thus eliminated,, data accumulation could be resumed. By 
choosing angles between 0.0 and 10.0 degrees, various portions of the diode 
array could be brought into the path of the beam. Proton flux was measured 
using an Ortec surface barrier detector (SBD) mounted a few centimeters in 
front of the faceplate. Sheets of Polaroid film were exposed to obtain beam 
area and flux density (particles/cm^-sec) , 



11 


Although Intended to operate with 512 fully Independent diode channeUi 
the Digicon could be tested effectively with a small fraction of these 
diodes active. Ten of the 512 diodes available were monitored through the 
use of a charge- sensitive hybrid preamp connected to each diode through a 
multipin header at the base of the tube. Each premnp was followed by a 
chain of electronics contaitving a pulse shaper, amplifier, discrimator, rate 
limiter and 16-bit count accumulator. Pulse-height analysis from any diode 
channel was possible, using switchable analog channels placed just after the 
preamps. These were followed by spectroscopy amplifiers and several multi- 
channel analyzers. As many as three diodes were monitored simultaneously, 
allowing measurements of coincidence rates between any pair or all three. 

It should be noted that although < he primary beam energies were 85 and 
40 MeV, passage of the protons through a 10-mil tantalum window used to 
isolate the beam line, as well as 0.5 m of air, the surface barrier detec- 
tor, and a few low-densicy light shields, reduced the energies of the pro- 
tons entering the window to 80.7 and 34.2 MeV, respectively. Passage 
through the 3 mm MgF2 window further reduced this to 75.4 and 21.5 >teV. All 
results will be referred to by these last values. 

Calibration 

Shown in Figure 5 is a pulse-height distribution of photoelectrons 

collected by a typical diode and produced by exposing the MgF 2 /CsTe face- 
plate to a UV penlight located several centimeters away. The peak of the 
distribution falls around channel 90, with the FWHM around 30 channels. 
Although each photoelectron has a kinetic energy of 22 keV, some spreading 
of the collected charge is expected on the basis of surface backscattering, 
and photoelectrons striking diode boundaries, both of which result in less 
Chan the full energy being deposited inside the diode. All pulse-height 
distributions presented here have Che horizontal axis calibrated with 
respect to Che phoCoelectron peak, which is assigned the value of 0.972 fC 
(6077 elect roir- hole pairs), assiming 1 electron-hole pair/3.62 eV of energy 
deposited in Che silicon, 


Pulse Height Analysis 


The differential pulse height distributions of signals obtained from a 
typical diode during direct proton bombardment of the Digicon diode array 
are indicated by curve A) in Figures 6 and 7» for 21«S and 75.4 MeV, respec’ 
tively. We have plotted the number of counts N(Q)dQ versus Q, where Q 
is the charge collected at a diode, Incident flux was low (10*^ - 10^ 
p+/sec-cii? , measured in real time with a Ortec surface barrier detector) and 
beam energy was sufficient (>8 MeV) to allow the protons to penetrate and 
completely pass through the 400 M silicon wafer with little energy loss (1 
MeV). No radiation damage was expected at these low fluences. Ihe shape of 
the differential pulse height distribution is typical of the response seen 
in all the diodes, and has been observed previously by us during proton 
bombardment of a similar diode array at other energies (ref. 19) and by 
others using low energy electrons and CCD's (ref. 10). 

A more useful representation can be obtained by plotting these data in 
integral form, indicated by curve B) in Figures 6 and 7. Here we show the 
total number of events in which the collected charge was greater than Q, 
versus Q. Two immediate conclusions can be drawn from these spectra. 

Using the 21.5 MeV result as an example, we see that the largest value of 
charge occuring in the spectrum (11 £C) is more than an order of magnitude 
larger than the amount of ionization that would be produced by a 21.5 MeV 
proton passing through the 2.5 M depletion layer. This result suggests that 
a substantial amount of charge released in the substrate beneath a diode 
reaches the diode depletion zone. A similar effect has been observed using 
low energy electron bombardment of thinned CCD arrays (ref. 10) , and users 
of surface barrier detectors have observed pronounced diffusion collection 
from depths as great as 500 p when these devices were irradiated with 
electrons (ref. 25). 

Second, by recording the value of the integral spectrum at the noise 
threshold (0.5 fC), one obtains the total ntinber of protons detected by the 
diode for the duration of the experiment. We have found this number to be 
larger than the number of protons incident on the area of the diode by 
almost an order of magnitude (for 21.5 MeV). This suggests that the ioniz- 
ing particle does not have to directly strike a diode to produce a signal in 


13 


chat diode. The net effect if to inctease the effective area or croaa~fec~ 
tion of each diode to several times its actual dimensions (which is 8.5 x 
10"5 cn? for the Digicon diode). Both of these effects become less pro- 
nounced at higher proton energies as a consequence of the lower specific 
ionization produced in the bulk, 

The shape of the differential pulse height spectra in Figures 6 and 7 
substantiates these conclusions. Most of the counts appear to be associated 
with small charge values, presunably originating from particle tracks cre- 
ated exterior Co and far from the diode boundaries. Since this region en- 
compasses a greater area than that of the diode itself, the protons, which 
strike Che array at random, have a greater probability of striking outside 
the diode, resulting in the observed distribution of counts. The largest 
charge values, occurring above the broad peak, presumably represent charge 
collected from particle tracks penetrating near the center of the diode, a 
relatively infrequent event. However, the factors contributing to Che exact 
shape of the pulse height spectra are complex, and will be discussed in 
detail in ’’Analysis." 


Coincidence Measurements 

The close proximity of one diode to another, as well as the observation 
that Che effective areas of adjacent diodes overlap, suggests that diffusing 
charge can leak to several diodes within a very short time interval, and 
thus create signals in adjacent diode pairs nearly coincidentally. This has 
indeed been observed in a 250-U thick CCD array exposed to highly ionizing 
cosmic rays, whereby dozens of pixels were simultaneously activated (ref. 
11). To study this effect, signal lines were connected from various combi- 
nations of diode pairs into a coincidence circuit with a one microsecond 
acceptance window. This unit produced one count for every current pulse 
which occurred in both diodes within one microsecond of each other. It was 
expected chat the coincidence rate would decrease with increasing diode 
separation, as the diffusing charge would have to disperse over a larger 
area, resulting in less charge reaching both diodes. 


The coincidence response of one diode with respect to another has been 
defined by dividing the total number of coincidence counts occurring between 
the two diodes during an arbitrary time interval by the average total number 
of counts obtained from both diodes during the same interval. Shown in 
Figure 8 is a plot of this fraction as a function of diode separation, for 
both proton energies used. Also shown are the results for 21.5 MeV and a 40 
microsecond acceptance window. As expected, the coincidence response falls 
of with increasing didode separation; the largest values always occurring 
for adjacent diodes, although there is a measurable effect up to six diode 
spacings apart (for 21.5 MeV protons). Also indicated in Figure 8 is the 
strong energy dependence of the coincidence. This is due to the specific 
ionization loss of the proton being less at higher energies, thus producing 
less charge per unit length to be shared among the diodes. At energies 
below 21.5 MeV, larger coincidence fractions are to be expected. 

Since the coincidence response for a given pair of diodes is presumably 
proportional to the extent to which each diode's effective sensitive area 

overlaps the other's, it should be possible to establish a dependence of I' 

this response on the geometrical arrangements of the diodes on the array. | 

To experimentally verify this dependence we performed a two- parameter 

analysis using the IBM 360/44 data acquisition system available to users of j 

the University of Maryland cyclotron. In this experiment, signals from a | 

• I ^ 

pair of diodes were each sent to separate analog- to-digital converters. If 
the conversion in one AOC occurred within 40 microseconds of that occurring 
in the other, both events were registered as a coincidence and sent to a 
dual parameter and display controller where the events were plotted and then 
stored on magnetic tape. A Z- intensified coutour display of several diode 
pairs is shown in Figure 9. 

In each display, the amount of charge collected in one diode is plotted 
versus that collected in the other, and only those events vAkich occurred in 
coincidence are plotted. The intensity of the data points is proportional 
to the number of counts each point represents. The figures suggest that 
adjacent diodes display somewhat different coincidence patterns Chan diodes 
separated by one diode. Based on our interpretation, the pattern observed 
in Figure 9 is not in any way fundamental, but depends on geometrical fac- 
tors such as diode shape, relative spacing of each diode on Che array, etc. 
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ConsequenCly , a full discussion of the nature of the coincidence pattern 
must be postponed until "Analysis," when we use the diffusion model to 
discuss the results presented here* 


* 

i 


DIFFUSION MODEL 
Introduction 

As described earlier, part of the kinetic energy associated with a 
penetrating charged particle is absorbed in the n-Si substrate of the diode 
array, creating hole*electron pairs. Ihese ionization- generated charge 
pairs then diffuse away from their point of generation, some of wliich reach 
the depletion regions of the reverse-biased diodes. Ihis chapter examines 
the physical theory behind the noise sensitivity of the diode array target 
as determined by hole diffusion and the discrete nature of the diode array. 

. Charge Creation 

We consider protons entering the front surface of a silicon diode array 
at normal incidence. This surface is divided into two regions: exposed 

diodes and the Si-Si02 interface (see Figure 1) . An energetic proton pene- 
trating the face of the silicon wafer will generate in the substrate, 
through collective plasmon excitations and subsequent thermalization, elec- 
tron-hole pairs at a rate of 1 e-h pair/3.62 eV of kinetic energy lost by 
Che proton.’ Initially confined Co the path of the proton, Che charge pairs 
thermalize within a few picoseconds and spread to a radius of about 0.1 p. 
The charges Chen immediately begin Co diffuse in a direction opposite Co the 
concentration gradient and disperse throughout the substrate of the wafer. 
For those which reach the edge of a diode depletion zone, the holes are 
swept into the depletion zone and the electrons are ohmically relaxed 
through Che substrate contact. The fact of electron hole pairs that reach 

the Si-SiO^ interface depends on the boundary conditions at this surface. 
Occasionally, a proton passes through both a diode's depletion region and 
the substrate beneath, where charge produced in both regions add quickly, 
producing a rather large current pulse; Protons penetrating the wafer out- 
side the boundary of a diode enter only the substrate resulting in less 
charge reaching a diode. An exact solution of the diffusion equation, to be 
given shortly, shows the collection currents peak very rapidly after 
generation in the bulk and the collection process is complete within micro- 
seconds. 
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Other minor factors affecting the amount of charge collected include 
the electronic time constants of the input circuit to the diode, as well as 
dielectric relaxation effects in the equivalent circuit of the semiconductor 
aubstrate-'-depletion zone. We shall ignore these effects here for clarity 
of presentation, cs they are significant only for large charge collection 
times. 

Finally, mention should be made of a phenomenon discovered in the past 
year, known as "field-funneling." Recently, Hsieh et al., (ref. 26) pres- 
ented definitive evidence that when an alpha-particle penetrates a p~n junc- 
tion, the generated carriers drastically distort the junction field. After 

the alpha-particle penetration, the field. Which was originally limited to 

* 

the depletion region, extends far down into the bulk silicon along the 

length of the alpha-particle track and funnels a large number of carriers 

into the struck junction. After a few nanoseconds, the field recovers to 

its position in the normal depletion layer, and, if the track is long 

enough, a residue of carriers is left to be transported by diffusion. A 

first order analysis by Hu (ref. 2,7) shows the funneling depth is equal to 1 

+ times the depletion region width for an n+/p junction (p-type sub 

strates), where M and p are the drift mobilities for electrons and holes, 

respectively, in silicon. For n-type substrates, the ratio becomes 

p /p , This result does not depend on charge density per cm along the track 
P ^ 

and should apply to proton tracks as well. 

Boundary Value Problem 

We begin by deriving an expression for the amount of charge reaching a 
diode from the substrate due to the passage of a charged particle. Making 
use of a Cartesian coordinate system, we approximate the silicon wafer by a 
homogeneous medium of infinite extent in x and y, and of finite thick- 
ness. In this and subsequent discussion, we shall assume this region to be 
of n-type conductivity. The front surface of the wafer, (z * 0), is par- 
tially implanted with diffused planar junctions, each of which is assumed to 
be a perfectly absorbing surface of zero thickness and finite area, (The 
details of the junction fabrication process do not concern us here). The 
remainder of the surface consists of the semiconductor-oxide interface 
boundary (Si-Si 02 ), where we assume incomplete charge absorption takes place 
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due to the effect of fast aurface atatea. Incomplete abaorption alao exiata 
at the back aurface of the wafer (e ■ W), although thla boundary nay be 
prepared under different conditiona than the front aurface> 

Aa aeen in Figure 10 t the origin of a rectangular coordinate ayatan ia 
centered on the front aurface of a typical diode with a rectangular croaa** 
aection of arbitrary dimenaiona. The -z axia pointi into the ailicon bulk» 
Faaaage of a ionizing particle, auch aa a proton, perpendicularly through 
the ailicon wafer at t ■ 0 createa, nearly inatantaneoualy i a line of charge 
parallel to the -z axia. The denaity of charge geiierated along thia track 
ia generated non-uniform and will be formally repreaented by the function 
p(z^). The coordinates of the charged particle's point of entry ia (x^, 

0), that of any point along the trajectory is (x^, y^, ). We aaaune 

normal incidence for simplicity; experimentally this is easily realizable. 

The volume density p of minority carriers (holes) in the quasi-neutral 
semiconductor medium under low injection conditions is governed by the time 
dependent diffusion equation given below (in C&rtaaian coordinates) 


' 3 p / 32 p 32 p 

— - + D 1 H 

,3t P\3 j^ 3j^ 



p(z ) 6(t) 
o 


( 1 ) 


where p (z^) Is the source function, the diffusion coefficient for 

holes in n-Si, and Xp is the hole lifetime. This equation gives the 
minority carrier concentration at the location r ■ (x, y, z) for a line of 
charge distributed along a track at r' “ (xq, y^) and created at t ■ 0. 

In general, the solution to this equation must satisfy several boundary 
conditiona at both the front and back surfaces of the wafer. 

For the front surface: 




z"0 


z«0 


V p 

o*^n 


V p 

in 


(at Si- 02 ) (2) 

(at depletion layer) (3) 
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For th« back surfaces: 


3p 

- - D ~ 
z'H " 9t 


z"W 


V p 


(4) 


where Vn is che surface reconbinatlon velocity at the respective bounds* 
ries. 


An analytical evaluation of the diffusion process in a mosaic target is 
quite involved. For a silicon wafer of finite thickness W, (1) becomes a 
boundary value problem with little chance of an exact solution. The reason 
for this la the need to specify three different homogeneous boundary condi- 
tions (of the third kind), one each for the diode depletion zone, the sur- 
face of the silicon wafer (z - 0) and its backside (z ■ W) (see Appendix 1). 
In fact, an exact solution would require detailed knowledge of the shape of 
the depletion regions and their geometrical arrangement on the array, 
Problems such as these with mixed boundary conditions (i.e,, Dirchiet and 
Neumann boundary conditions specified on the same surface) are notoriously 
protracted but have been attempted by a few authors (refs. 28, 29, 30), As 
von Roos points out (refs, 28, 29), solutions are usually in the form of a 
series expansion of orthogonal functions appropriate to the coordinate 
system used, where the expansion coefficients are obtained by solving a set 
of dual inhomogeneous integral equations. To the author's knowledge, a 
solution to the particular problem posed here has not appeared in the dif- 
fusion theory literature. 


Approximations 

Mixed boundary values are necessary in theoretically evaluating experi- 
ments where surface effects predominate due to the close proximity of the 
charge source to the surface, such as the EBIC work described by Von Roos 
(ref, 28), However, simpler boundary conditions may not suffice even when 
the charge source lies deeper in the medium. For the case of long particle 
tracks penetrating deep into the silicon wafer, most of the charge is 
created away from the surface influences. However, it remains to be seen 
whether or not surface states continue to affect the amount of charge 
reaching the diode depletion zone. Nevertheless, there appear to be a 
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ceet«in con<Utioni under which it leemi reasonable to relax boundary condi- 
tion (2), and aaaume (3) applies over the entire surface. This approxima- 
tion should adequately describe devices having a large fraction of their 
surfcce covered with diodes (such as CCD arrays). It should also hold for 
linear diode arrays and short particle tracks if the particle track is 
created well within the diode's boundary, away from the influence of surface 
effects at the nilicon-oxide boundary. Also, as recombination velocities 
are usually quite high at depletion zones, (>10^ m/s), due to the intense 
electric fields present at the junction, it is customary to let «, 
which is equivalent to demanding the charge density vanish at z ■ 0. Under 
these two approximations we are thus assuming the entire front surface to be 
a perfectly absorbing surface. 


Solution 

As shown in Appendix I for a medium unbounded in x and y, and infi- 
nite or semi- infinite in the z-coordinate, the Greens function solution to 
(1), under the approximations just covered, has the form 

G. “ — i — exp /- -f l ... - ~V / exp (-n^ D t} K(n,z) K(n, z^)dn (5) 

P 4irD^t \ ADt T y a P ® 

P ' P P' 

where cZ - (x - Xo)2 + (y - yQ)Z 
and a“0 0<z<» (setu- infinite) 

a«-» -«<z<«» (infinite) 

and the K(n, z) are eigenfunctions for the homogeneous problem correspond- 
ing to the boundary condition at z ■ 0. 

For a median of finite thicknebs the integral is replaced by a sunma- 
tion over eigenfunctions appropriate to the boundary conditions at z •• 0 and 
z - W. 

This solution was obtained by solving (1) using an impulse source 
function for a point along the trade 
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( 6 ) 


G ■ g 6(x-x ) 6 (y-y ) 6 (z-z ) 6(t) 
s 0 0 o 

where g is the source strength at that point, and whose value we have 
chosen to be unity for the time being. Later we will integrate Gp over 
the charge density function p(Z(,) to obtain the final result. 

By means of an example, the kernel K(n,z) appropriate for an infinite 
medium, K(n,z) ■ exp(inz) when substituted in equation (5), gives 


G ■ f(r,t) / -1 exp {-r? D t - in(z-z )}dn 


-• 2s 


resulting in the Greens function for an infinite medium 

+ (z-z ) 
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P (4sD t)3^2 
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exp - 


4D t 
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t 

T 


(7) 


For a semi- inf inite medium* in which V * at z ■ 0 the appropriate 
kernel is K(n,z) ■ 2ir sin(nz), when substituted in equation (5) gives 


(4s D. 


i e«p/- 

[ 




Cz-z )^' 

exp I 2^1 *expl- 


4D t 
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/ (z+z )2\ 
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I 
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f \ A r A / . 


( 8 ) 


The result, although obtained formally by integral transform techniques, is 
physically motivated by use of the method of images whereby the semi- infi- 
nite medium and a charge sink located at + Zq. Iherefore, the charge 
density vanishes at z ■ 0, which is equivalent to demanding an infinite 
recombination velocity there. 

Introducing a boundary at z *• W requires a different kernel to be 
substituted into equation (2). In this case, K(n,z) is given by 


K(n, z) 


$2 + (v /D )2 
n w p 


, w(B^ + (v /D )2) + V /D 
^ n w p ' w p_ 


1/2 


sin 0 z 
n 
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for V ■ • at z ■ 0. This results in the following; expression 


G 

P 


■ exp (- I exp D t} 

Air Dpt \ 4Dpt Tp/n-1 

02 + (v /D )2 

X 2 _Ji — E sin 0 z sin 0 z 

w(0n + /D )2) + V /D ■ " " ° 

^ n w p • P 


where are the positive roots of 0 * cot(0W) + Vy/Dp ■ 0. 

For the limiting cases of + 0 and Vy • at z ■ W this Greens 
function reduces to 


AirD t 
P 


exp - 


V 4D t T / n«l 

\ n n/ 


exp I- 


(n - 1/2)2 ^2 0 

P 


w • / X . / Iv ir 

X sin (n - «.) — z sin (n ^ —) — z 


2 W 


2 W 


( 10 ) 


for Vy > 0 (perfect charge- reflecting back surface) with a similar expres- 
sion for (perfect charge- absorbing back surface) with (n -■ •^) tar 

placed by n. 

The current density at a point (x, y, 0) on the surface is given by 


j(r,z^,t) 



( 11 ) 


Next, the current flowing into the entire diode from the point source 
^*o» yo> "2o) is given by integrating equation (11) over x and 
y, where the integration limits corresponds to the physical dimensions of 
the diode. 
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Upon integrating this expression over the tracklength L, we obtain 
the total current entering the diode resulting from charge diffusing away 
from the entire trajectory 

L 

I(r,t) ■ / p(Zq) i(r,z^,t) dz^ (13) 


where p(zq) is the generated ele Ttron-hole density per unit length. 

Since p(zq^ usually non-uniform we may choose for the purposes 
of numerical computation to approximate the integral in equation (13) by a 
summation over thin segments or layers 


I(r,t) 


N 

I 


n*l 


/ i(r,z ,t) dz 
o o 

z 

n 


(14) 


where zj^ ■ 0, z^ ■ L and N is the number of segments used to approximate 
P(zq) . We assume Pj, is uniform over each segment having thickness 
8 ■ *n+l ~ *n‘ average value of pjj is given by 


p * e * (dE/dx) /w 
n n 


where 

e ■ electronic change 

w ■ energy required to create 1 e-h pair 
(dE/dx)n ■ stopping power for nth layer 

Finally j the total charge Qg collected by the diode from the sub- 
strate track is obtained by integrating over all time 


00 

Qg - / I(r,t)dt 
o 




dydxdz dt 
0 

z“0 


(15a) 
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The integrations indicated in (I5a) may be carried out analytically for 
both semi- infinite and finite solutions, with only the time integral requir- 
ing numerical integration. This results in the following expressions for 

Qe 


N 

Q. ■ I p„ w 

* n-1 " " 


(15 b) 


t. 

where 



for a semi- infinite medium, and 
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— f X£ (erf) exp 
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- cos 
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W 


(17) 


for a wafer of thickness W having a finite recombination velocity v^ 
at z ■ W, and where 
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and where 


■ 4D t, a„ ^ 3 W, b ■ V W and L . ■ D t 
p ’ m in ’ w dp p 

where ot^ are the positive roots of oi*cot(a) + b ■ 0 

Written this way, Wq represents an "effective track length" or the 

length of a track 'seen' by the diode with lOOZ efficiency, as opposed to 

the much smaller efficiency with which the actual track is seen, which we 

shall call the "charge collection efficiency" C , where C « W (z .,-z 

^ n n n n+l n 

We may write (I5b) in terms of Cq as 
N 

Q„ ■ I q_ C (18) 

n"i 

where qq is the charge created along the segment of track having length 
Zn+l “ *n* 

To equation (15) we must add Qq, the charge collected from the 
depletion layer, the surface diffused layer, and the "funneling zone." Qq 
may be written as 


where 


■ charge density/unit length along W^ 

P 

W . (1 + -^) (2epM V ^)^^2 ( 19 ) 

o n o 

n 

and the physical parameters are defined in Table 1. Note Qq “ 0 if the 
proton enters the wafer at a point outside the perimeter of the diode. 


Semi-Infinite vs. Finite Approximations 

In the semi-infinite approximation, the total charge collected from 
the substrate depends only on the diffusion length , plus geometrical 
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Table 1. Physical characteristics of digicon diode array. 


SUBSTRATE 

Doping 

Resistivity (p) 

Thickness (W) 

Dielectric constant (e) 
Hole Mobility (Pp) 
Electron mobility (m,^) 
Diffusion coefficient (Dp) 
Bias voltage (Vq) 

DIODES 

Type 

Area 

Pitch 

Width 

Length 

Depletion width 
Funneling depth 
No. diodes on array 


n - type 
4.5 ohm-cm 
400 M 
12-13 

70 cm^/sec-v 

1380 - 1400 cm2/sec-v 

14.8 cm^/sec 

+5 V 


P'fn diffused (0.7 )-' 
with Boron 

8.5 X 10“® cm^/diode 
50 p 

40 p 
250 p 

2.5 p 
1 P 
512 


factors involving diode dimensions and track length, Ihe effect of the 
water back surface on the flow of the excess charge carriers -is ignored. As 
a general rule, if the wafer thickness W >> Lj this approximation is valid 
to a high degree, Another important geometrical ratio is that between diode 
size and W, which if small enough provides an independent criterion for 
ignoring the back surface. The importance of the magnitude of these geomet- 
rical factors cannot be overstated, for although the functional dependence 
on them appears complex, dimensionless ratios among them may be formed which 
shows more clearly the minimum number of variables upon which charge collec- 
tion depends. 

For a Digicon-type diode and substrate, both expressions equations (16) 
and (17) are plotted as a function of Ld and shown in Figure 11 for a charge 
track directly beneath the diode (x^, - ■ o). Wte have normalized the 

values plotted to the track length W, which has been set equal to 400 p, 
corresponding to the physical situation in which the proton penetrates the 
wafer and thus track length and wafer thickness are the same. Notice we are 
assuming p(a^) Is uniform along the entire track. As might be expected, 
increasing the diffusion length beyond the track length results in little 
additional charge collection, as evidenced by the leveling off of the curve 
around Lj » 300 U . These curves are practically identical in both the 
serai- inf inite and finite approximations, for both values of v„, except 
for Ljj w, where the curve for v^ » “ is slightly depressed. 

Actually, only charge originating from that portion of the track near- 
est the back surface is affected by conditions there, as can be seen- in 
Figure 12 where we show the collection efficiency from four 100 p track 
segments plotted versus Lj for Vy » 0 and Vy ■ As the curves show, 

very little charge reaches the diode from the last 100 P anyway, since the 
collection efficiencies barely exceeds 0.5%. 

Although the amount of charge collected by the diode from a given par- 
ticle track depends strongly on the diffusion length and the track depends 
strongly on the diffusion length and the track length, the maximum charge 
collectable from a track is ultimately determined by the diode area. This 
effect can be seen clearly in Figure 13 where the effective track length is 
plotted versus actual track length, for a number of values for L<j. Notice 
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that beyond a limit defined by the envelope of all the curvea, additional 
increases in either Lj or W fail to further the amount of charge 
collected by the diode. In particular, the maximum theoretical charge 
collection efficiency from a 400 M track is about 9 %. 

Summarizing, we find that unless Lj is of the same order as W, 
back-* surface effects can be safely ignored and the much simpler semi-infi- 
nite approximation used in modeling charge diffusion from particle tracks. 
Second, if the wafer thickness is much greater than front surface diode 
dimensions, the semi- infinite model is valid irrespective of Lj, since 
the diode subtends vanishingly small solid angles front points far along the 
track. For the case of the Digicon diode/substrate, both considerations 
work to our advantages; consequently we shall employ the semi-infinite 
approximation in our modeling. 

Diffusion Current Waveforms 

Although we are primarily interested in the total charge collected at 
the diode, it is of some interest to examine the diffusion current waveforms 
in order to explore the physical soundness of the results of our derivation. 
We confine our attention to the Digicon type diode- substrate whose relevant 
physical parameters are sunmarized in Table 1. Ihe general shape of the 
diffusion current versus time waveforms is illustrated in Figures 14, 15, 16 
for a track created by a proton entering the wafer at y^ “0, “ 0» 50 and 

100 y , respectively (using the same coordinate system defined at the begin- 
ning of the chapter and in Figure 10). These curves are plotted using 
equation (13) and the semi- inf inite Greens function with a track length W 
of 400 m , as before. We do not show drift currents due to charge collection 
in the depletion layer or the funneling zone, as these occur on subnanosec- 
ond time scales. Evident is the strong dependence of the shape and magni- 
tude of the current pulse on the distance ftom the proton's entry point to 
the center of the diode. In addition to the expected decrease in current 
with increasing distance of entry, both a peak shift and width increase is 
apparent. The time delay arises naturally as a result of the longer dis- 
tances over which charge must diffuse to reach the diode. This also 
produces in a smearing out of the arrival times of each charge carrier, due 
to the random nature of the diffusion process, resulting in larger peak 
widths. 
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Decreasing the diffusion length also has a pronounced effect on 
both the shape and magnitude of the current, as indicated in the eeriea of 
figures beginning with Figure 17. Since Lj controls the depth to which 
the diode "sees" into the silicon bulk, it is to be expected that smaller 
values of L^j will lead to a reduction in the magnitude of the current, 
as well as the time spread of the current pulse. This latter effect follows 
from the smaller time differential for charge arriving at the diode from the 
beginning and end of the track, resulting from a decrease in the effective 
length of the track, (see Figure 13). 

Since the magnitude of each current pulse is extremely small, (nc more 
than 100 pA at the peak for a 20 MeV proton), and lasts for less than a 
microsecond, it is unlikely that diffusion current waveforms due to particle 
tracks can be measured with much precision at present. Fortunately, much of 
the relevant information can be derived from the integrated charge. Experi- 
ments obtaining this information were described and, in the next section, 
those results are interpreted in light of the model presented here. 
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ANALYSIS 


Introduction 

Our primnry task in this chapter is to account for the experimental 
frequency distribution of noise charge observed during irradiation of the 
Digicon diode array with 21.5 and 75.4 MeV protons, and presented in "Exper** 
imental Work." These pulse-height spectra of collected charge represent the 
accumulation of tens of thousands of separate events each consisting of a 
proton entering the wafer at some random point (xq, y^) and generating 
charge along a straight line path. The electron-hole pairs subsequently 
spread in all directions, inducing a transient diffusion current in a p-n 
junction diode at the surface of the wafer. Thus, all pulse-height spectra 
simulation algorithms in this area are faced with two fundamental physics 
problems; charge creation and charge transport. 

Once charge is created at a point in the wafer, the subsequent trans- 
port and collection of that charge at a diode is a rather complex diffusion 
theory boundary value problem but one which can be solved under the simpli- 
fying approximations outlined in "Diffusion Madel." However, it is well 
known that the ionization losses of fast charged particles fluctuate, even 
if the particles are monoenergetic, due to the statistical nature of the 
energy” loss process. Thus, the number of charge carriers creat>%d in the 
silicon by each proton is not constant, but characterized by a probability 
function we shall refer to as the energy- loss straggling distribution. This 
means the magnitude of charge collected by a diode from each event is deter- 
mined not only by diffusion, but also by statistical aonsiderations involv- 
ing the ionization process itself. 

To clarify this point, consider the magnitude of the charge Q at the 

endpoint of the integral pulse height spectra in "Experimental Work." This 

point, Q , represents th'e largest amount of charge collected by the 
max 

diode for that particular energy. For the 21.5 MeV data shown in Figure 6, 

Q is around 11.2 fC, which is equalivenC to about 70,000 charge carriers 
max 

collected. At 75.4 MeV this point occurs at around 5.4 fC . From symmetry 
considerations, these events should be due to protons penetrating near the 
center of the diode's collecting area, where the largest solid angle is 
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available to the ionization track. Diffusion theory, however, predicts a 
maximum collectable charge Q of only about 7.1 fC from a 21,5 MeV proton 
track, (assuming a uniform p(zq) and no statistical fluctuations from 
ttacl; to track). Furthermore, as will be illustrated shortly, calculated 
pulse height distributions using this assumption result in spectra which 
poorly approximate most of the upper half of the experimental distribution. 
Clearly a more rigorous approach is needed. The next section outlines these 
ideas in more detail. 

Diffusion-Predicted Pulse*Height Spectra 

Before we introduce energy loss fluctuations, let us examine the degree 
to which diffusion theory alone accounts for the experimental results. To 
do this we musi: provide some familiarity with the simulation techniques used 
in this work and the assumptions behind them. To simplify matters, we begin 
by plotting as a function of x^ and y^ the two dimensional charge col- 
lection efficiency matrix C„ defined in expression (18), for a uniform 
track of length W * 400 u. This is shown in a three-dimensional representa- 
tion in Figure 20 for a Digicon type diode/substrate, and Lj ■> 90 p. As 
expected, the largest values corresponds to tracks located beneath the 
center of the diode, falling off most nearly exponentially with increasing 
distance from the diode boundaries. (The peaked structure in Figure 20 has 
no special significance and is due to the particular shape of the collection 
area of the Digicon diode). When each value in this matrix is multiplied by 
the total charge generated in the silicon by a proton (not including the 
depletion zone), Figure 20 represents the amount of charge collected by the 
diode from the substrate, plotted as a function of the proton's entry point. 
A computer program can then be used to sort these values into a frequency 
distribution which represents a rather crude approximation to the experi- 
mental distribution of collected charge. 

This calculation was carried out and is shown plotted in Figure 21 in 
integral form alongside the experimental result for 21.5 MeV protons, while 
Figure 22 shows the 75.4 MeV result. In both figures curve (a) corresponds 
to Wq « 0 (no depletion layer) , while (b) corresponds to Wq “ 3.5 M 
(given by expression (19) in "Diffusion Model"). In both cases the lack of 
a statistical treatment results in a gross under stimat ion of the maximum 
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charge collected, and faila to account for the overall ahape of the diftrl** 
bution for large values of charge. Fairly good agreement is aeen in ther 
lower half of the spectrum, at both energies, which is also that portion of 
the spectrum whose shape is independent of Wq, This suggests diffusion 
theory alone suffices in accounting for charge collection from distancti 
particle tracks, or at least tracks not directly underneath the diode col- 
lection area. 


Energy Loss Fluctuations 

A better approximation can be constructed by allowing for fluctuations 
in the amount of ionization produced in the silicon from event to event 
(track to track). In the previous treatment we simply multiplied each value 
in the collection efficiency matrix by the mean energy loss. A, given 
by the familiar expression due to Bethe [see for example, Rossi, (ref. 44)] 


0.30058 s 
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wher.'e Z, A, and I are, respectively, the atomic number, atomic weight, 

. C , 

and mean excitation energy of the medium, 2 ~ is the shell correction, and 

Z 

5 is the density effect correction. 


where 
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on mass where me? is the electron rest energy and 0c the velocity of the 

particle. For small 0, e is also very small. Because of the small 

ndx 

value of ra/M for protons or mesons, the factor in square brackets is 

practically unity, except when 0 1. For 0 1, + — , which means 

that almost the entire energy of the particle can then be transferred in a 


single collision, (We have neglected the shell correction and density ef- 
fect correction in (20) as they are insignificant for our range of 0). 
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Nov the energy loit by a charged particle in a thin abaorber is subject 
to fluctuations as previously mentioned. Landau (ref. 31) derived a trans~ 
port equation for the straggling function for the energy losses d small 
compared to the initial energy of the incident particle: 


0 

" /** w(c) f(A - e,s)de - f(A,a) w(e)de (21) 

9s o o 


(b ■ A for A < e j b ■ e for A > e ) , 
max’ max max 


where f(A,s)dA denotes the probability that a proton in penetrating a 
distance s in the target, will experience an energy loss between A and 
A -f dA as a result of collisions with atomic electrons; and w(e)de is the 
differential collision cross-section for single collisions with an energy 
loss e. 

Thus, the mean energy loss A is given by 

A ■ /* f(A,s)A dA 

o 

A rigorous solution to equation (21) has been given by Vavilov (ref. 32), 
who showed that the shape of f(A,s) was characterized by the ratio of A 
to Gmax* denote by K, where 

X - 0.30058 -- s/e - -i— . (22) 

02 A e 

max 


If each of the energy- transfer to electrons is small compared to A, 
which is the case if K » 1, f(A,s) is approximately Gaussian, given by 

f(A,s) exp {-(A - A)2 k/2%^ (I - 02/2)} (23) 

5 — (1 - 02 / 2 ) 

(see also Jackson (ref. 33)). 

If the total number of atomic collisions is small, which occurs for 
both thin absorbers or fast particles, (l.e., x « 1), then individual 
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coUilionf can make significant contributions to and the distribution 
is asymmetric, with a tail corresponding to large values of This is the 
so-called Landau distribution, given by 


lim 

K a- 0 f (A , s) 


i- /** e" cos [x(log X + X)J dx, 
ir5 0 2 


(24) 


where 

X ■ X /k - log K »■ t — £ - 1 - + Y - log tc, 

^ 5 

Y " 0>577216 (Euler's constant). The shape of the Landau distribution 
is indicated in Figure 23. 

For intermediate values of k, the rather complex but exact solution 
due to Vavilov applies, and is given in Appendix 2. Both the Landau and 
Gaussian distributions are limiting cases of the exact Vavilov solution. 

For the Digicon diode/substrate all of these types of fluctuations contrib- 
ute to the overall shape of the p'ulsc^* height spectra in Figures 6 and 7. In 
addition, at the lower proton energy (21.5 MeV), this shape function is 
dominated by fluctuations in the 3.5 |i depletion 1 aye r/funne ling zone, where 
K «< 1 and the above theoretical treatments are not valid. 


Computer Program PULHGT 

To correctly simulate the observed pulse height distributions, the 
above considerations force us to treat each incoming proton as a separate 
event. A computer program called PULHGT was developed in order to acconr 
plish this task. PULHGT performs the calculations indicated in (18), and 
using Monte Carlo techniques to determine x^, y^, and q^, calculates 
pulse height distribution of charge which is collected as the result of N 
protons striking the diode array at random. This is accomplished by first 
calculating the threct-dimensional charge collection efficiency matrix 
C^(x^,y^, z^) for a large number of points representing the volime in the 


35 


substrate to which the diode is sensitive. This matrix is stored in PULHOT 
and during the simulation each point in this matrix coirespondir)|j to a 
particle track at x^, is multiplied by the charge created along each 
segment of the track having length - z^. The number of charge carriers 
created by each proton in either the depletion zone or along a segment is 
determined by computing' the energy loss straggling function f(A,s) appro^ 
priate for the region through which the proton is passings and using a 
random number generated on an interval from 0 to 1 to select values from the 
normalized integral of f(A,s). In this way, is calculated inde- 

pendently for each segment z„ in expression (18) as well as for the 
depletion layer and the process is literated for as many times as is 

necessary to build up to the equivalent statistical accuracy of the experi- 
mental distribution. 

In some situations it is convenient to allow the ionization to vary 
from event to event, but to retain the ass unction that the ionization 
density remains constant over the track length. In this case our simulation 
Is essentially two-dimensional, as we are calculating for each point x^, 

Yo an average charge collection efficiency for the entire 4Q0 M track, 
and multiplying this value by the charge produced along the track. (This 
was the procedure used in the diffusion-only simulations previously, where 
we neglected energy- loss fluctuations completely). 

This latter procedure was followed in building up the theoretical 
pulse height distributions shown in Figure 24. The Vavilov parameter for a 
400 u track is large enough (>1) to assume a Gaussian distribution of energy 
losses in the substrate. For the 3.5 M depletion layer/funneling zone, ic ■ 
0.0297> meaning we should use the Landau distribution to determine energy 
loss in this region. The program PULHGT computed the energy of the 40 MeV 
protons after passing through the 10 mil tantalum beam port and the 3 mm 
MgF 2 window also using a Gaussian distribution of energy losses. Two random 
mxnbers were used to determine Xq and yo for each proton. 

We show the corvesponding predicted and observed integral/differential 
distributions at 21.5 MeV in Figure 24. A trial-and-error procedure was 
used to determine the value of L^j needed to produce the best fit. This 
turned out to be Lj * 90 p . We have normalized the theoretical integral 
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dlstrlbuClon to the experimental dlatrlbutlon at a point around 3# 5 fC» as 
this seems to produce the best overall £lt< From Figure 24 ve see that 
agreement between theory and experiment Is fair except for the very low 
«1»4 £C) charge events, In which case the simulation underestimates the 
amount of charge collected. Notice also that most of the large charge 
events (>6 fC) appear to be due to energy-loss fluctuations In the depletion 
layer. 

When a similar simulation is performed for 75.4 MeV protons, the agree- 
ment between the integral distributions is excellent until around 1.2 fC, 
where the two curves dramatically part company (see Figure 25). Unlike the 
previous result at 21. 5 MeV, here the fit Is poor for both Wq * 0 and Wq 
■ 3. 5 M . Evidently, energy loss fluctuations in the depletion layer do not 
fully account for large Q events at this higher energy. 

Three-Dimensional Simulations 

Explaining the discrepancy in the 75.4 MeV integral spectra at large 
charge values requires rejecting the assumption of a uniform p(zq). 

Tnitlally we assumed a uniform density of electron hole pairs along the 
particle track to simplify the derivations and expedite the calculations, as 
full three-dimensional simulations require substantial computing time. The 
reason this assumption must be rejected lies in realizing that the diode 
"sees” different segments of nearby particle tracks with different degrees 
of efficiency, determined mainly by two factors: the geometrical solid 

angle subtended by the diode at a point on the track segment, and of course 
the diffusion length. (By 'nearby' tracks we mean those occurring directly 
beneath the diode collecting area.) Thus, the diode is somewhat sensitive 
to large fluctuations in the charge density along these tracks, and this in 
turn affects the number of large charge events seen. Of course, these fluc- 
tuations are more severe at higher proton energies, since the number of 
collisions per unit length is smaller. We show the dependency of on 

Zq for a track centered at xq * yo =0 in Figure 26. 

Charge collection from distant tracks is less influenced by variations 
l>^p(zo). Indicated in Figure 27 is Cn(zo) for a track entering at 
*0 ■ 100 M , y ■ 0. Peak sensitivity occurs around Zg ■ 50 y , dropping 
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rather sharply on both sides. Only the portion of the track centered on the 
peak is seen by the diode, compared to the larger fraction seen for nearby 
tracks, This analysis suggests why diffusion theory alone suffices in 
determining the experimental charge distributions for low to median Q (due 
to distance tracks), but fails in the high Q part of spectrum (due to 
nearby tracks). 

To simulate the 75.4 MeV data we chose to divide the track into 10 
segments each having a length of 20 M for a total track length of 200 M. 
(Although the actual track length is 400 p , we have shown that for L^] < 

150 p , very little charge is calculated to reach the diode from points along 
the track beyond 200 p). This choice was a reasonable compromise between a 
desire to adequately sample the varying collection efficiency along the 
track, and budget the expensive computing time. The theoretical spectrum 
obtained for 75.4 MeV protons is presented in Figure 28. Although a slight 
deviation remains for very small Q, the overall agreement between theory 
and experiment is much improved. 

At 21.5 MeV the fit is not as precise (Figure 29). There occurs a 
discrepancy for very small Q which we have attributed to surface effects 
not adequately modeled here. The rest of the calculated spectrum follows 
its experimental counterpart about as well using the three-dimensional 
treatment as it did with the two-dimensional approximation. This further 
indicates the insensitivity at low proton energies to fluctuations along the 
substrate track, in striking contrast to the 75.4 MeV result, in which sub- 
strate fluctuations dominate the spectrum. 

Low Q: Surface Effects 

A possible reason for the improved fit at 75.4 MeV for events occurring 
below 1.4 fC may be seen by considering where the protons which cause these 
events enter the diode array. Using the program PULHGT we plotted the 
values of Xq and yo responsible for causing counts within a 'window' 
of charge values occurring in the pulse height spectrum, whose width and 
possible could be varied at will. Such a plot is shown superimposed on the 
diode array pattern in Figures 30 and 31 for 21.5 and 75.4 MeV protons, 
respectively. These plots show where protons had to enter the wafer in 
order to result in charge between 0.5 and 1.4 fC being collected by the 
center diode. 
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At 21.5 MeV most of the events originate in particle tracks created 
well outside OlOO p) the boundary of the center or 'reference' diode. This 
is a region in which much of the charge roust diffuse under the Si/Si 02 
interface portion of the diode array pattern in order to reach the diode. 
Much of this charge will reach the surface before getting to the diode 
depletion zone, and according to the boundary condition Vq ■ be 
completely absorbed there. If this boundary condition is not valid then 
some charge would be 'reflected' at this interface and reach the diode. 

This would probably result in more counts being observed at low charge 
values, but we do not attempt the calculation here. (This was partly our 
reason for not normalizing the theoretical spectrum to the experimental one 
at the noise threshold at 0.5 fC) . 

However, we remarked in "Experimental Work" that if protons enter the 
array within the 'protective' region near the diode elements, where surface 
effects should be minimal in determining the amount of charge reaching each 
diode, then our assumption of a completely absorbing surface at z ■ 0 may be 
a valid approximation. This is borne out in Figure 31, since most of the 
protons? at 75.4 MeV causing counts in the above region of interest apparent- 
ly enter the array close to the center diode's boundaries. 

High Q: incorrect Distribution Function? 

Returning to the problem of large charge values in the spectrum, we 
have found for charge values above 3.2 fC (at 21.5 MeV) that aside from 
being confined to the interior of the diode, there is little correlation 
between the location of the "charge window" and Xq, y^. Actually, the 
region near the endpoint of the spectrum represents relatively infrequent 
events when large energy transfers occurred in the depletion zone and sub- 
strate. This effectively increases the charge collection efficiency beyond 
that predicted by expression (16), explaining the discrepancy described at 
the beginning of the chapter. 

Unfortunately, theoretical treatments of energy loss in extremely thin 
absorbers, such as the depletion layer, are not yet consistent with experi- 
mental results. Vavilov/Landau theory does not include the effect of 
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resonant collisions with distant atomic electrons, which can significantly 
alter the shape of f(A,s) in these regions (ref. 42). Our use of the Landau 
distribution in this region is therefore probably incorrect, and may explain 
the poorer fit seen for large Q at 21.5 MeV. Corrections to the Vavilov 
theory, most notable that given by Blunck and Leisegang (ref. 34), have been 
shown not to agree with experiment (ref. 38), although recent efforts appear 
to be more successful (ref. 39). Vavilov's theory works very well in the 
range of K used at 75.4 MeV to model the substrate p(zg) (ref. 40, 41). 

We assume this is responsible for the excellent fit seen there. 

Analysis of Coincidence Results 

An important aspect of charge diffusion from particle tracks created in 
multi-diode thin silicon chips is the near simultaneous triggering of sever- 
al diodes from a single proton event. Every diode in the vicinity of a 
proton track sees the track with differing degrees of effectiveness, and at 
slightly different moments in time, depending on the solid angles and diffu- 
sion time constants of the situation. This effect, besides providing an 
opportunity for further testing the validiaty of various physical models, 
can have many practical consequences with respect to the use of thse devices 
as paraticle detectors. Some proton-induced noise-rejection schemes to be 
employed when the Space Telescope Digicon diode array is flown in orbit 
depend on various anti-coincidence decisions whose effectiveness is deter- 
mined by the magnitude of the coincidence effect. 

We have not been as successful in modeling the coincidence and effec- 
tive area results indicated in "Experimental Work," as these phenomena 
depend rather strongly on the number of low Q events collected, vdiich our 
model underestimates at low proton energies. We do see qualitative agree- 
ment, however, indicating the physical principles involved are probably 
sound. In the limit of infinite front surface charge absorption, we predict 
crosstalk and enhanced collection is reduced by 50% over the absorbed 
amount. Our results suggest, at least, that these effects may be controlled 
by altering the surface preparation of the wafer during manufacture. 
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Summary and Conclusions 

The combined problem of generation/ transport of ionizing-radiation 

created Charge carriers in a silicon microelectronic device has been treated 

here for the first time. Given the extreme complexity of this process, 

conventional wisdom would dictate that a useful simulation procedure is 

practically impossible by any except Monte Carlo methods. Indeed, Monte 

Carlo techniques have proven successful in modeling energy loss distribution 

(ref. 35) as well as charge transport simulations in integrated circuits 

(refs. 13, 14). However, a full scale Monte Carlo simulation of the 

combined generation/transport problem would probably tax the limits of most 

computer installations. From this viewpoint, it is encouraging that the 

« 

predominately analytical approach pursued in this paper, with the 
unavoidable simplifying assumptions, has met with reasonable success. 

Furthermore, we have clearly delineated the areas in which the model 
fails to work, by establishing a connection between those areas and the 
assumptions which underly them. This is most important if Monte Carlo 
modeling is attempted in future work, as this approach requires a clear 
understanding of the physical processes occurring in the system, and their 
relative importance. Two important problem areas brought out in this work 
are surface effects and energy loss in thin depletion layers, the latter 
being of critical importance in simulating proton-induced charge noise in 
digital circuits having very low 

Finally, we believe the interaction of protons with these devices 
presents a challenge to both the experimentalist and the theoretician by 
providing a simple tool with which to study charge creation and transport in 
a semi-classical environment. Future devices with their even greater levels 
of integration will surely require quantum mechanical modeling, and the 
considerations raised here will continue to be relevant, if not dominating. 


APPENDIX I 


GREEN FUNCTION SOLUTION TO DIFFUSION EQUATION 


In the absence of electric fields ^ the continuity and current density 
equations for holes in n-type silicon are: 
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where p„ ” number density of holes in holes/cm**3 
Sp “ hole generation rate 
Up * hole recombination rate 
Jp ■ hole current density 
Dp “ hole diffusion constant 
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Under low injection conditions (i.e., where the injected carrier densi- 
ty is much less than the equilibrium majority carrier density), Up can be 

approximated by (p - P„„)/t , where p is the thermal equilibrium minor- 

n no p ^ no ^ 

ity carrier density and Xp the hole lifetime. 

For the three-dimensional case under low injection conditions, equa- 
tions (I.l) and (l.2) combine and reduce to 
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For an instantaneous line source of charge created at t * 0, we have 
Sp * p(zQ)6(t), where p(zo^ generated charge per unit length. 
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For t > 0 (1.3) is homogeneous and we may substitute p ■ (p - p^^) exp(t/Tp) 
leaving, in Cartesian coordinates 


P \3)^ 3y2 3 z 2 / 3t 


p(z^) 6(t) 
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The Green function G(x,x') appropriate to (I. 4a) satisfies the in 
homogeneous equation 


-D + L-£. + L-2.\ + i£. ■ 6(x-x ) 6(y-y ) 6(z-z ) 6(t) (1. 4b) 

P\a,^ 3y2 at o '.o o 

Equations (1.4) are familiar from the theory of heat conduction. 

Ozisik (ref. 43) has developed Green function solutions to equation (1.4) 
under a nunA>er of different boundary conditions, using Fourier integral 
transform techniques. We may define the Fourier transforms needed to solve 
equation (1.4) as follows. For a medium infinite in x and y and 
infinite, semi- infinite or finite in z, the transform is defined as 


pO ,V ,n ,t) ^ K(n,z) p(x,y,z,t) dxdydz (1.5) 

!■<)> III^D 

along with the inverse transform for a semi- inf inite or infinite medium 


p(x,y,z,t) - /* /" / e K(n,z) p(0,v,n,fc) dBdvdn (1. 6a) 
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where 0 semi- infinite 0 < z < * 

o ■ 

-• infinite < z < ® 

and for a finite medium 0 < z < W 


o» 00 “ "i(Px + ^y) 

p(x,y,z,t) ■/ / I K(n__,z) p(0,v,n_,t) d3dV (l.6b) 
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The function K(Tim,z) is called the transform kernel and the Tim's 
are the eigenvalues. The type of kernels and the equations for the eigen- 
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values depend on the combinations of boundary conditions at z " 0 and z ■■ 
and are tabulated on page 50-52 and 481-483 in Ozisik, 

Solution 

We first remove the time dependence in (I, 4b) by applying the Laplace 
transform to both sides of the equation, giving 
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P 3x^ 3y^ 3z^ 


- s ■ 6(x-x ) fi(yy_) 6(z-z ) (1.7) 
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where 


G » L(G] and L[6(t)] - I 

Applying the transform defined in equation (1.5) to both sides of 
equation (1.7) we get 
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a simple algebraic equation whose solution is 


i(0x^ + vy^) 

-e K(n,z^) 

D (02 + v2 + p2) + g 

P 


( 1 . 8 ) 


Applying the inverse Laplace transform to equation (1.8) yields 
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K(n.p ) 
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Finally, applying the remaining inverse transform (I. 6a) to (1.9) and 
using the integral formula 
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we obtain for a semi- infinite medium 
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and Gp ■ g exp(-t/xp). 
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APPENDIX II 


VAVILOV SOLUTION TO LANDAU TRANSPORT EQUATION 


Details on Vavilov's treatment of energy loss fluctuations can be found 
in the literature (ref. 32). His solution to the transport equation is as 
follows: 


f(A,s) dA -1 ♦ (X .K.$2)dX 

5 




jc ^K(l ^ 02y) ^ ^Kf, ^ 


X - - K(1 + S2 . y) 

^ e 

max 

Y ■ 0,577216 (Euler's constant) 
f|(y) ■ S2[iog y + Ci(y)] _ cos y _ y Si(y) 


f2(y) ■ ytlog y + Ci(y)I _ cos y _ y Si(y) 


Si(y) • ^ du (sine integral) 

o u 

Ci(y) ■ /y ^ .S 5 . . H du (cosine integral). 

— flO y 

Tabulations of the Vavilov distribution have been performed by Seltzer 
and Berger (ref. 36) for protons and mesons, and the computer code used in 
their numerical integration was graciously supplied to us by Steve Seltzer. 
Modifications of this code were necessary in order to calculate the integral 
of f(A,s), which is performed in the program VAVINT. The original code 
may be found in NBS Report 9498 (ref. 37). 
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APPENDIX III 
COMPUTER PROGRAMS 

These are the FORTRAN codes used in the simulations presented in this 
paper. DIODE generates the charge collection efficiency matrix C,^(Xo, 

Xo* given the circuit geometry and substrate diffusion length. The 

integration indicated in (16) uses IMSL routine DCADRE, which performs nu" 
merical integration of a function using cautious adaptive Romberg extrapola- 
tion. The integrand in (16) is contained in the function routine XYZTL. 

The Cn matrix is written onto a disk file for use by the simulation 
routine PULHGT. 

Simulations are carried out by transporting a proton through the wafer, 
layer by layer, and computing the charge generated in each laye:, Charge 
reaching the diode is determined by multiplying this dE/dx vector by the 
charge collection matrix C^ as in (18), Proton entry point (x^, y^) 
is generated using random numbers, and the energy loss in each layer is 
determined by sampling the appropriate probability density function f(A,s), 
The latter is accomplished in the following way, First the density function 
is converted into a integral distribution function P(A,s) given by 

P(A,s) « I - /^ f(A',s)dA (111,1) 

—00 

If f(A,s) was properly normalized, P(A,s) will be uniformly dis- 
tributed on the interval (0, 1), and be a proper distribution frop which 
to sample f(A,s), Actually, this can be accomplished analytically only if 
(IIl,l) can be solved for A in terms of P(A,s), Since this is generally 
not possible, numerous interpolation schemes are employed to sample from 
tabulated values of P(A,s), These values are generated in program VAV, 

I 

where both f(A,s) and P(A,s) are calculated for either Vavilov or 
Blunck-Leisegang theories. 

The actual simulations are carried out in PULGHT, which contains suf- 
ficient documentation for understanding. These programs were written on a 
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5 

I 

I 

f 

PEC-IO uBxng version 5 of FORTRAN 10, end version 8 of the Internstional 
Mathematical Subroutine Library (INSL). The graphics program PLTGEN was 
written from PLOT- 10/ ADVANCED GRAPHING II, Release 1,2, and from the LEVEL I 
version of the 4010A01 PLOT 10 TERMINAL CONTROL SYSTEM, Release 3.3, avail- 
able from TEKTRONIX. 

The interpolation subroutine INTERP was taken from Bevington (ref, 45), 


I 
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CALLING, PHOCRAK tO CKEATE CHAWn COLLECTION EFflCIENCV 

MATRIX/ <iSED IN PULSE-HEIGHT SI*<CLATIOhS 

* 

DIMENSION AyT0TC3‘^/121/lf )/DI0Dn(3,4) 

COM PON /P Aha M / X .> / y 'i , DL, T H I C K/ A/ B / R et a 
COMMON /LlMnS/XMllvX.MAXy yhIN/XMA X 
EXTERNAL XXZTL 
DOUBLE PRSCISION FtU«E 

DICICON DIODE DIMENSIONS (MICRONS) 

DATA ( (DIODE (I/J),lsl/3)/iJ=l/4)/-2G,-2'’/»6*858/2U/2f-/ 
X6.860/-7b,i25/75,75,i7b/12S/ 

CALL nEfiSST(''/Li:VOLD) I SUfRESS IMSL ERRCR PRINTOUT 

TYPE 1C*A 

PaRPAT(XX/lX'KNTBR DIFFUSION LENGTH/ X INC/ V INC;*$) 

ACCEPT */ DL/ IC/ ID 

CALL CNFIG(XyTOT/DIODE/IC/ID/SIEP) 

TYPE 30? 

FORKATdX/lX'SNTER OUTPUT FILE NAMBt'S) 

ACCEPT 2.?/ FNAME 
FORMAT (A6) 

OPE N( UMT=21 / FILE =FNAm£) 

WKITE(21/») STEP 

*RITE(21/2jfO(((XYTOT(l/J/N)/I=l/35)/J=l/X21)/A«l/ir.A 
FOR MAT (3 SEX 0.4) 

CLOSE! UNIT=21/FILS=FKAP.E) 

CALL EXIT 

END .. 


I 
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C 

4 

3 

2 

1 


creates CHAPGr CGLL5.CTI0N 2Ff’ICI2«CV ><ATRIX Cn<XO/YO/Zn) 


SUbROUTIfiE CNF IG 1 ('X¥ TOT /DIODE /IC/ ID/ STEP) 

DIWEWSION DICDi(3/4) 

DlhSNSIOts X.YTOT(35/121/l;) 

COMPaN /PARA.M/XC/Yr/DL/THICK/A/H/BiiTA 
COkFON /LlhlTS/XRIN/XFAX/VFlN/YFAX 
EXTERNAL XYZTL 
PI =3.1 41 59265 4 ; AA=^ * ; U =4DC 

CONsl./SClRT(16*PI) 1 semi-infinite NORMALIZATION 

a=^:.;b=2C. 


THICK=«-A I 
XO=0.0;Y'*=-12v. 

STEP=3. I 
DO 1 J=1/121/ID I 
DO 2 1=1,35/ IC I 
DO 3 N=l/10 I 
DO 4 M=l/3 I 
XMIN=DI0DE(M/1) 

XMAX-DILjBE C M,2) 
YMIN=DI0DE(M/3) 
YMAX=DI0DE(M/4) 


SEGMl-NT THICAMiSS 

RESOLUTION (MICRONS) 

INCREMENT YO 
INCREMENT X?) 

INCREMENT ZC 

DIGICON DIODE HAS THREE PASTS 


NUMERICAL INISGRATION OVER TIME FROM AA -> 0 USING DCADRE, 
IMSL ROMUERG INTLGRATION ROUTINE 

XYTCT(I,J/N)=XyTOTa/J/N) + DCADRE(XYZTL^AA/0/1.0E-3/C./E/liN) 
CONTINUE 

A=A THICK; Ei:H 4 THICK 
XYTOT(I/J/{0=XVT(JT(I,J/N)*CON/THICK 
CONTINUE 
A=C .‘''; 3=1 HICK 
XL=xrtSTSP»IC 
• CONTINUE 
XD=C.O 

¥0 = V**fST£P*ID 

CONTINUE 

RETURN 

END 
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FUNCTION XYZTL(P) 

SSMI-INFINITE CHARGE COLLECTION EFFICIfcflCY 
integrand for NUHrJRICAL INlEGRATION OVER TIME 
£RF(X) IS AN IMSL CALLABLE FUNCTION 

CONNON /FARAN/)(^i/YO/DL/ThlC«/A/B/B£TA 

COMMON ZLIMITS/XMIN/XMAX/YmIN/YMAX 

ARG1=(XMAX-Xn)/R 

ARG2s(XMIN-X0)/R 

ARS3=:( YMAX-Y0)/R 

ARG4=(yMIN-Y0)/P. 

EHF1=ERF(ARG1) 

ERF2=ERFCARG2) 

eRf3=SRF(ARG3) 

ERF4=ERF(ARG4) 

X=£HF1-ERF2 

Y=FRF3-ERF4 

2=LXP(-(A/R>*^2) - EXP(-(B/R)**2) 
T=EXP(-(k/(2*DL))**2) 

XYZTL=X*Y*Z*T 

RETURN 

END 
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PROGRAM PfJLHGT 

P£RF0RmS simulatijks of pulse-hright data 

DIMENSION CN (35, 121/10 ),C(6),COUNTa50w) 

DIMENSION X‘(lSC0),DUAL(l,S/5) 

DIMENSION C4u(4)/C240(4),C85(4),AC(4),ATC<b),C14D(6) 
DIMENSION XC1(2),C840(2) 

dimension XLARL1(59)/XLABL2(S9)/ALAB’3(59)/YLAEL1(59) 
DIMENSION YLABL2(59) 

COMNON/ThlN/GP03(7D),GFO325(7O)/GPO35(?O) 
COMMON/THICK/GP17(40)/GP19(40 ),GP21 (40 ),GP23 (40 ),GP2 5 (40 ) 
COMVON/SiMPL£,/GTHlH(10C)/GTHICKaOO>/NOMl/MUN2 
CQMMQN/YVALUE/VL1(100)/YL2(100) 

COMMON /THINGS/ENRG/DSEEO/KMAX/NBQL/LLOA/LLDB/IiRSUMl/NSDH 
l/XA/VA/n,T2/LA/LE,NID,LLDl/LLD2/LLD3/LLD4/T/NSDK3/NhD/STEP 
2/NND1/NND2/NSTAT 

COMMON /PLOT/NPTS/XMlN/XMAX,YMIN/yHAX/Xl/X2/X3/X4/NP/KG/LT/ 
lIS7P/ISY«,ITX,m/IG 

DOUBLE Pk’ECISION FINAME/DSESD 
REAL NCFT'J ’ 

COIKCIDENCE DATA FROM DIGICQH TEST RUN 

DATA €40/(5.28/0.23/0.13,0.03/ 

DATA C240/0. 73/0. 61/0. 34/0. 22/ 

DATA C85/0. 20/0,077/-!. ,.014/ 

DATA C840/0. 36/11.29/ 

DATA XCl/1.,2./ 

DATA XC/1./2./3./6./ 

DATA XTC/1./ 2./3./4./5./6./ 

LABELS FUR TEKTRONIX PLOTTING 


ni /'D'/'e*/ 'r 


a^r/g/e/ /b#a/r/r/i/e/r^s/ /v/n^ 


DATA XLABLl/'N'/'u'/' 

• '/ 'r', 
b'o'/'l'/ 

l*e^/ 's*/ ')'/' ■*/*C'/*o*/*l'/*l'/'e*/'c*/'t*/'e 
DATA yLABLi/'C'/'o'/'u-/-n'/'t'/ 's',53*' V 


'/'d'/16*' '/ 


DATA XLAHL2/ 

I'e', 'c', 

'C 

m 

•/ 

• «K • * a ‘ 

/ n / a 
'd'/43** 

/ r / 

'/ 

g / 

9 ^ # # 4 

e / 

9 9 

i ^ / 

0 / 

'1', 

'1' 

DATA XLABL3/ 
l'r',47*' V 

'D 

'/'i'/'o' 

'/'d'/' 

6 / 

' '/'N' 

p 4 

# U 4f 

'p V 

'b',‘ 


. DATA YLABL2/ 

'C 

'/'O'/'l' 

'/'n'/' 

c / 

'i'/'d‘ 

/By 

'n*/‘ 

'c'/‘ 

•e* 


1* ^/^F*/ *r*/*3*/*c^/*t*/*i*/*o*/*n*/39’ 


V 


HEAD IN VAVILOV INTEGRAL PROBABILITY DISTRIBUTIONS 
CALL READIN 

MRITE(5/*) (yLl(I)/GTHlN(I)/I=l/N0Ml) 

TYPE 40 

ACCEPT 41/FINAfcE 
OP£N(UNIT=23,FILE=FINA«S) ' 

READ IN charge COLLECTION EPFICIENCk MATRIX Cn 
READ (23/*) step 

HEAD (23/100 /END-70) (( (CN (1/ J, KA/I^l, 35), J=l, 121 ), Kil/lO) 
TO CLOSE(UNIT=23/FILS=FINAM2) 
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xoe FORMAT CaSElC.^) 

4C FDRmATC input fils WAM2? ') 

41 FDRMAT(AiC) 

LAyEK,=K-l 
TYPE 97 

97 FORVATdX/lX'ENTER SPECTRUM, RESOLUTICN (E-H PAIRS):'S) 

ACCEPT NHOL 

71 TYPE 5C1 

501 FORMATCIX/IX'ENTER PROTON ENERGY < NINCOW (MEV): 'S) 

ACCEPT */ ENERGY 
TYPE 93 

93 FORWATdX/lX'RSJCCT PROTONS NOT ENTERING AREA CEFINEC 

1 By -X/X,-y,y:'s) 

ACCEPT XAMII,/XAMAX,YAMIN,Y AMAX 

TYPE 43 

43 FORMATdX/lX'FCR PRIMARY SPECTRUM ENTER LLf)#ULC:'$) 

C 

C LLD: LOWER LEVEL DISCRIMINATOR/ OLD: UPPER LEVEL DISCRIM. 

C 

ACCEPT */ LLDA/LLDb ‘ 

TYPE 44 

44 FORMATdX/lX'FCR PSCTOh ENTRY POINT MAP ENTER LLD/ULD;'$) 
ACCEPT *, LLD1/LLD2 

TYPE 46 

46 FORMATdX/lX'FOR COINCIDSNCS DSTECTICN ENTER DIODE NC. 

1/LLD/ULD:'$) 

ACCEPT */ HWD/NWD1/NWD2 
TYPE 92 

92 FORMATdX/lX'STOP SIMULATION WhEN TOTAL COUNTS 

1 BETWEEN XI AVD X2 EXCEED Y:'S) 
accept */ LLD3/LLD4/NIT 
TYPE 69/ LAYER 

69 FORMATdX/lX'ENTSR THICKNESS OF DEPLETION LAYER AND EACH OP 

1 '12' LAYERS; '5) 

ACCEPT */ T1/T2 
TYPE 89 

89 FORMATdX/lX'ENTER WHICH LAYER AND HOW MANY DICDES;'S) 

ACCEPT »/ LB/NID 
LA=L0 

IF(LB.EQ.O) LAsl 
IF(LB.EQ.C) LB-LAYER 
TYPE 888 

888 FORMATdX/lX' NAME FOR PROTON ENTRY MAP FILE;'$) 

ACCEPT 41/ FINAME 
OPEN(UNIT=21/FIL.E=FINAmS) 

WRITE(21/*> LLD1/LLD2/NIT 
TYPE 910 

910 PORMATdX/lX'STRAGGLING ?:'$) 

ACCEPT */ NSTAT 
TYPE 887 

887 FORMATdX/lX'TYPS 1 FOP DEBUG') 

ACCEPT *, T 
C 

C INITIALIZE PARAMETERS FOR RINDOm MATERIAL AND 

C DETERMINE ENERGY LOSS OF PROTON PASSING THRQUGfc DIGICOfI WINDOW 

C 

ZA=.454;ZI = .09088;DEWSE=2.6C'l I LITHIUM FLOUFIDC 

BET A=BTA(ENERGY) 

CALL IUN1Z(RETA/30CC.,2A/ZI/DENSE/EWC/DWQ,£MAX/DUM/DUM) 

G- 
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C ZERO accumulator ARRAYS 

C 

DO 1 K=1/1SC0 
COUNT(K)=C 
1 CONTINUE 

DO 2 K=l/5 . 

DO 2 J=l/5 

2 DUAL(l/K,J)-0.0 
DO 3 J=l/6 

3 C(J)-0.0 
C 

C BEGIN PROTON IRRADIATION SIMULATION 

C 

DSEEDS121212.D0 
AMiNXsXAMIN - XAHAX 
AHINYsYAMlN - VAMAX 
DO 33 HPROTsl/lOOOOOOO 
XAsGGUBFS(DSEED)'*fAMINX + XAMAX 

<AsGGUfiFS(DSSED)*AMlNY + YAHAX 

C 

C, DETERMINE ENERGY OF PROTON AFTER RINDOte BUT BEFORE DIODE ARRAY 
C 

DELh=GAUSS(EK0,D40,DSE£D) 

ENRG=ENERGy-(DELW/100D.)t ENERGY AFTER WINDON 
C 

C INCREMENT PULSE HEIGHT SPECTRUM 

C 

CALL PULHGT(XY^C, COUNT/DUAL) 

IF(NS0M3.EU.NII) GO TO 659 I JUMP OUT OF ITERATION 

33 CONTINUE 

659 CLOSE(UNiT=21j-FILS=FINAME) 

C 

C SUMMARIZE PARAMETERS AND RESULTS 

C 

A3-ENERGY-(EHO/1000) 

2=14;A=28.085;ZA=Z/A;Z1=.17C?;DENSE=2.328 1 SILICON 

BETA=BTA(A3) 

CALL LANDAU (£ETA,T1/ZA/ZI, DSN SE/DElr£l#ULAM,V,Y,V) 

CALL LANDAU(BETA/T2,ZA/ZI/DENS£,EE2/E2yDLAM,V/«/V) 
WTHK=3.G;DWT=Tl;bT=T2 

AI = ENERGY;A2=WTHK;A4sDA0/10C)0 ;A5sDNT;A6=E1 
A7-E1*1000/3.62;A8=D£1;A9sDE1*1000/3.62;a 1&=BT 
All=E2;Al2=E2*100y/3.62;A13sDE2;A14sDE2»10f>0/3.6 
A15S4.72E2* NFP.CT;Al6=FLaAT(NSUM)/FLGAT(NRSUMl) 

C 

C DEFINE OUTPUT ABOVE. ..TYPE OUT BELON ' 

C 

TYPE 999/Al/A2/A3,A4/A5,A6/A7/A8/ A9/A10/A11,A12,A13, A14/A15,A16 
999 FQRMATCIX/IX'INCIDENT PROTON ENERGY: 'F4.1' MEY' 

1/lX'AFTBR 'F2' MF NlNDON: 'F5.1' MEV +/- 'F6. 4' MFV'/ 

2' AVERAGE ENERGY LOSS 'IX/IX'IH 'F5.1* MICRON DEPLETION WIDTH: ' 
3F6.2' KEV 'FS" S-H PAIRS'lX/3t7X '+/- *F7.3' KEV 'F5'E-H PAIRS'lX 
4/' IN 'F5.1' MICRON BDLK:'F7,.2' KEV 'F7* E-H PAIRS'1X/17X'+/- 
5 'F7.3' KEV *F6' E-tl PAIRS' 

6/' TOTAL DOSS: '3' P+/CM**2'1X/1X- 
7'EFFECTIVE AREA/DIQDS AREA; 'F5.2) 

TYPE 883/ LLDA/NSUM/LLD3/NSUM3/NRSUM1/NPROT 
883 FORMATdX/lX'TCTAL COUNTS ABOVE '15' RERE ' I7/1X/ 

IIX/'TGTAL COUNTS ABOVE '15' WERE '17/lX/lX/ 


1 ASSIGN SEED 10 GGUEFS (DSEED) 


I ITERATION LOOP 
i GENERATE RANDOM X 
1 GENERATE RANCOH Y 
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2TOTAL PROTONS TO HIT DIODE ItERE ‘’I5/1A/1X 
3'MUMBSR OF ITERATIONS kAS *16) 

TYPE 88A/XAMIK^XAMAX/YAPIN/YAMAX/NUMU#I 
884 f'OR^ATaX/lX'NUPDEP. Cf PROTONS TO HlT kINOOW DEFINED 6V * 
14FS' WAS *16) 

C 

C PLOT PULSE HEIGHT DISTRIBUTIONS BELOW 

C 

Xl=l.;X2=l.;X3=7.76;X4s5«75;NG='yES';NP=*YSS* 

LT=-4;ISTP=1 

XSYN=0;ITX=l;ITY=2;IG=2;NPTS=KHAX-l 
TYPE 10010 

lecio FORKATCIX/IX'GRAPH PULSE HEIGHT CISTRIBOTI-K ? 1=YES') 
ACCEPT GRAF 

1F(GRAF.SQ.0) 60 TO lOCll 
DO 8 Isl^KMAX 

8 XCDsFLOATCNHOL^d-D) 

45 TYPE 50 

50 FORFATdX/lX'SNTER XFIN,XMAX,YMIN/YMAX-) 

ACCEPT XMIM,XrtAX/VMIN/YMAX 

CALL PLTGSN(X/COUNT) 

CALL LAaL(1.9,0.25/'H',XLABLl,0.vl*5/l. 5/10^0) 

CALL LABL(0. 2/3. ,*’V'/YLAbLl/0. /I. 5/1^5/10, P) 

CALL TINPUKJK) 

CALL ERASE 
TYPE 121 

121 FORMATdX/lX'GRAPH AGAIN? 1=YBS') 

ACCEPT */ GRAG 
IF(GRAG.EQ.l) GO TO 45 
IF(FLAG.SQ.l) GO TO lOCll 
C 

C OUTPUT PULSE HEIGHT SPECTRUM TO DSK FILS 

C 

TYPE 40 

accept 41/ FINAME 
□PEN(UMIT=:21/FILE=FINAME) 

WRITE(21/*) (XCK)/COUNT(K)/K=l/NPTS) 
CL0SE(UNIT=21/FILE=F1NAMS) 

C 

C PLOT INTEGRAL CF SPECTRUM 

C 

TYPE 1122 

1122 FORHATdX/lX'PLQT INTEGRAL OF SPftCTRUM— T YPE IS'S) 

ACCEPT »/ INT 
IF(INT.NE.l) GO TO 10011 
FLAGsl;SUMrt=0.7 
DO 1124 I-NPTS/1/-1 
SUMM=SUMM ♦ COUNT(I) 

COUNTdlaSUMM 
1124 CONTINUE 

GO TO 45 
C 

C PLOT DUAL PARAMETER; X VRS. Y VRS. CCIKCIDSNCE INTENSITY 

C 

l&Cll TYPE 51 

51 FOKMATdX/lX' iiISH TO PLOT DUAL PARAMETER?*) 

ACCEPT */ GRAG 

IF(GRAG.SQ.O) GO TO 56 
ITYsl;X3s6.;X4-6. 

913 TYPE 52 
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ACCEPT HDIODS 
TYPE 5^ 

ACCEPT XMINyXKAX/VMIN/y«AX 
DO 1(! 1=1,5 
CnDNT(l)=(I-l)»10v0 
DO 1C J=l/S 
n=c 

IF(DUAL(KDIODE,I/J).GT.l) HsALOGlO (DUAL(HDI0DE/I, J)) * 0.5 

xa)=(j-i)*iooo 

isvmm;npts=i 

CALL PLTGEN(X/C0U?{T) 

NP='iJO';NG='NO' 

1C CONTINUE 

CALL LABL(2.6/0.25,'H',XLABL2,C./1.5/1.5,10/C) 

CALL LABL(U. 25 , 2 .S/'V*/XLABL 2 / 0 ./ 1 . 5 /l.S/lC/ 0 ) 

CALL TINPDT(KJ) 

CALL ERASE 
NP='YES'';NG=-YES' 

TYPE 121 
ACCEPT *, GRAG 
IFCGRAG.SG. 1) GO TO 913 
C 

C PLOT COINCIDENCE INTENSITY VRS. X 

C 

56 UO 57 J=l,lOOO 
COUNT(J)=O.C 

57 CONTINUE 

DO 58 J=l/5 
X(J)=(J-1)*1000 
Se CONTINUE 

TYPE 52 

52 FORPAT(lX/lX'ChOOS£ DICDE 1 OR DICDE 2 :'$) 

ACCEPT *, NDIODE 

IF(NDIOO£.E«.0) GO TO 11111 
TYPE 53 

53 FORMATClX/lX'ChOQSE CHANNEL NO. TO PLOTZ'S) 

ACCEPT tiCHAti 

TYPE 54 

54 FORRATUX/IX'ENTER 1 TO NORMALIZE Z '$) 

accept *, NMCRP 

DO 55 1=1,5 

COUNT(I)='DUAL(NDIOD£,NCHAN,I) 

IFINNORM.Eu.l) COUKTa)=DUAL(NDIODE,RChAN,I )/CCONT(l) 

55 CONTINUE 

NPTS=5;KP=‘’YES*;NG='VSS*;LT=0>ISTPal;ISVM=0 
TYPE 50 

ACCEPT *, XMIN,XMAX,YMIN,YM/X 

X3=7.75;X4=S.75 

TYPE 59 

59 FORFATdX/lX'LIN Y (1) OR LOG Y (2) ?Z'$) 

ACCEPT *, ITY 
CALL PLTGEN(X, COUNT) 

CALL LABLC3. 4,0. 25, 'U'’,XLABL2/C . ,1.5,1. 5,10, C ) 

CALL LA8L(C. 25, 3., 'V',YLABLl/0. ,1.5,1. 6,10,0) 

CALL TINPUT(P) 

CALL ERASE 
GO TO 56 
C 

C PLOT COINCIDENCE FRACTICNS VERi^US DICDE POSITICN 

G- ■ 
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11111 LT=-4;lTY=X>J!3=7.75;X4s5.75 

DO 151 11=1 /KID-1 
C(II)=C(II)/KSUM 

ISl CON TIN UC 

XMIN=l;X>'.AX=6;VKIH=0;YKAX=i. 

2323 NPTS=4;ISYM,=1;1STP=1;KP='YES';nGs*YES*; 

TYPE lOlM 

10101 FORFATdX/lX'LIN Y (1) OR LCO Y (2):'$) 

ACCEPT */ ITY 
lF(lTY.EU-,2) YHIN=.31 
CALL PLTGEK(XC/C40) 

NG=*N0';NP='N0*;ISYM=2 
CALL PLTGFH(XC,C247);ISYM=3 
CALL PLTGEN(XC/C85);ISYM=4;NPTS=N1D-1 
CALL PLT6EN(XTC/C);ISYW=5>NPTS=2 
CALL PLTGEN(XC1,C84C) 

CALL LARL(3. 6/0.25/ 'N'/XLAEL3/C./1.5/I.S/10/0) 
CALL LABL(0.25/2.3/'V'/YLABL2/0./l«5/1.5/lC/0) 
CALL TINPOT(JJ) 

CALL ERASE 
NPs'YES'ANGs'YSS' 

TYPE 121 
ACCEPT »/ GRAG 
IF(GRAG.EQ.I) GO TQ 2323 
CALL EXIT 
END 


SUBROUTINE PULHGT(XY/C/CQUKI/ DUAL ) 

THIS SUBROUTINE INCPSNENTS TriE PULSE-HtiGHT SPECTRUM AND 
PERI’ORmS coimcidfncs data SImULATICNS 

DIMEWSIDN CH (35/121/ 10), COUNT (1>/C(1) 
dimension DUAL(1,5/5)/DEDX(1- ) 

COMMON /THINGS/EWRG,DSEED,hMAX/KHQL/LLiiA/LLDfi/KRSUMl/NSUN 
l,XA/yA/Tl/T2/LA/LE/M10/LLDl/LLD2/LLD3/LLC4/T/NSUM3/NND/S'fEP 
2/NMD1/KMD2/NSTAT 

DOUBLE PRECISION DSEED 

FIND CHARGE DEPOSITED AT XA, YA / 2 

CALL CHARG£(DEPLT/UEOX,NSAV) 

COMPUTE COINCIDENCE AND DUAL PARAMETER DATA 


NX=0 

DO 299 K=1,NID 
DEPLET=0.0;SUM=0.v 
XANX=XA - NX 
yanx=-ya 

IF ((NX/50)/ 2*2. SU. NX/50) YANX = YA 
NX=NX SC> 

II=INT(AbS(XANX/STEP) + 1.5) 

JJ=INT(YANX/STEP ♦ 51.5) 

IF(Il.GT.35.0R.JJ.GT.121.0R.JU.LT.l) GO TO 299 
IF( K. EAI.NSAV) DEPLST=DEPLT 


57 


f 


ORIGINAL PAGE *8 


C 

C DETERMINE. CONTRIBOTXQN OF CHARGE FROM EACH LAVER 

C 

DO 23 LL=LA/L« 

SUMsSUM DEDX(LL)*Cftai/JJ/LL) 

23 CONTINUE . 

C 

C ADD DEPLETION ZONE CONTRlBOTIQh 

C 

TEST-SUM 4 DBPLET 

IF(K.E«.l) REFsTEST I SAVE VALUE OF DIODE 1 

lF(K.Ey.NWD) REFlsTSST I SAVE VALUE OF DIODE NWD 
I!f(K.Ey.l) GO TO 299 
C 

C ItICREMENT COINCIDENCE AND DUiiL PARAMETER COUNTERS 

C 

IF(EF;F.GE.LLCA.AND,TEST.GE.LLDA) C(K-1)=C(K-1) ♦ 1 
60 TO 299 

1D?REF/1000. + 1.S;.3D=1EST/100D. + 1*5 
IF(ID.GS.5) ID=S 
IF(JD.GE.5) JD=5 

IFtK.Eg.2) PUAL(K-1/ID/JD)=DUAL(K-1,1D,JD) + 1. 

299 CONTINUE 

C 

C REJECT IF MWD COUNTS OUTSIDE THIS tiINDOli 

C 

IF(REPl.LT.Nk01.0R.R£Fl.GT.Nl«D2) RETURN 
C 

C WRITE INFORMATION TO DISK FILE 

C 

IF(REF.GE.LLD1.AND.REF.LE.LLD2J hRlTE(21/*) REF/ 
1XA/VA,DSPLET/DEDX(LA) 

C 

C ACCUMULATE PULSE HEIGHT DISTRIBUTION 

C 

IF( FSF.GE.LLDA.AND.REF.LE.LLDP) KSUMsNSUN ♦ I 
IF(REF.6E.LLD3.AND.RSF.LE«LLD4> NSUM3sNSUH3 4 1 
K=REF/NH0L 4 1.5 
IF (K-KNAX) 4/4/2 

2 IF(K-1 500)15/3/3 

15 KMAX-K 

GO TO 4 

3 KslSCC’ 

KMA J?— K 

4 COUNT(K)sCaUNT(K) 4 1 

RETURN 

END 

C 

C 

C 

FUNCTION INS1DE(X/Y) 

C 

C THIS FUNCTION DEPENDS ON SHAPE OF DIODE COLLECTION AREA 

C HERE HE ARE USING UDT/ELECTROB VISION DIGICOK DIODE MASK 

C 

INSIDE=C- 

XABsABSCX);VAB=AbS(y) 

C 

C AM I INSIDE DIODE BOUNDAHIES ? 

(T 
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IF( XAB.LE.22 .5. AND. VAB.LE.77. S) INSIOEsl 
IF(XAB.LS.9.35.AND.y.GE.77.5.AND.y.L£.122.S) lASIDSsl 
IF(XAB.LE.22.5.AN0.y.GE.122.5.ANU.y.L£.I77.S) iNSIOEsl 
RETURN 
END 


SUBROUTINE CtiAHG2(DEPLT/DSDX/ NSAV ) 

THIS SUBROUTINE DETERMINES P«(Zn)— -CHARGE DENSIiy/SB GMENT 

dimension DSDXd) 

COMMON/THIN/GPC3(70)^GP'325(7C)/GP035(70; 
CONNQN/THICK/GH7(40)/CP19(4C)/GP21(40),GP23(4C),GP2 5(40) 
COMMON/SIMPLE/GTHIN(1CC)/GTH1CK(100),NUN1,NUM2 
COMMON/VVALUE/VL1(100),VL2<10D) 

COMMON /THIN6S/ENRG/DSEED,KNAX/NH0L,LLDA/LLDB,NRSUM1/NSUM 
l/XA,yA/Tl/T2/LA/LB/NID/LLDI/LLD2/LLD3,LLD4/T, NSUM3/NMD,STEP 
2^NND1/«V(D2,HSTAT 

DOUBLE PRECISION DSEEO 

INITIALIZE VARIABLES 


NTc.RMS=2/NPTS=75;NSAV=0 

INITIALIZE PARAMETERS FOR ABSORBING MEDIUM 
Z^=14>A=28.085;ZA-Z/A;ZI-.m;DENSEs2.328 I SILICON 

DID PROTON PASS THROUGH A DIODE? 

NX=0 

DO 5 K=1/NID 
XANX-XA - NX 
VAN Xs'" VA 

IF((NX/50)/2*2.EQ. NX/50) VANX^XA 
NX=NX 50 

I=INT(ABS(XANX/STSP) 1.5) 

J=INT(VANX/STEP -»• 51.5) 

IF(I.GT.35.0R.J.GT.121.aR.J.LT.l) GO TO S 
FLAG=INSIDS(XANX/YANX) 

IF(FLAG.EQ.l) GO TO 7 
CONTINUE 

IF NOT GO ON TO COMPUTE ENERGY LOSS IN BULK 
GO TO 337 

IF SO REMEMBER WHICH DIODE PROTON HIT 
j^SAV^K 

IFCK.EQ.i) NRSOMl-NRSUMl + 1 

IF SO COMPOTE ENERGY LESS IN OEPLETICN ZONE 

IFCTl.EQ.O) GO TO 337 

NPTS=NUM1 

.BETAssBTACSNRG) 

CALL I0NI2(BETA,Tl,ZA,2I/DENSE/E^DE/k,G^D0M) 
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335 

C 

C 

C 

337 


338 


10 


C 

C 


C 

C 

C 
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IF(NSTAT.EQ.n) GO TO 335 
Ul-GGUBFS(USSSD) 

CALL IWT£RP(GTbIN/YLl/KPTS/NTEHK&/Ul/TfeIK) 
FAClsG*W 

FAC2=1 + bSTA»»2 - C. 577216 ♦ AL0G(G> 

B=E FAC1*(THIN FAC2) 

DEPLt=S*10r/7./3.62 

CNixG'El^KG - £/lC’00. 


DETERMINE ENERGY LOSS IN EACH LAKER GF BULK 


NPTS=NUM2 
DO 10 K=1/LB 
BF.TA*BTA(ENP.G) 

CALL I0NIZ(BETA/T2/ZA/ZI/DENSE/E/DE/KVG,DUM) 
IF(HSTAT.eU«0) GO TO 338 
IF(C.GE'.l) E^GAUSSCE/DE/DSEED) 

IFLG.Gc;.!) GO TO 338 
UlsGGUBFS(DSEED) 

CALL INTERP( (.THICK/ VL2/NPTS/NTERRS/01/THIK) 

FAClsiG^W 

FAC2=1 + BSTA**2 - 0.577216 ♦ ALOG(G) 
e=E ♦ FAC1»(THIN -f FaC2) 

DFDnK>=En000./3. 62 
KT2=K*T2 

SNRGsENRG - E/IOOO. 

IF(T.EU.i) TYPE */ ENRG/XA/YA/KT2/DEFLT/DBDX<K)/THIH/G 

CONTINUE 

RETURN 

END 


FUNCTION GAUSS (5/DE/DSEED ) 

OUTPUT IS NORMALLY DISTRIBUTED AROUND £ NITH SlGHAsDE 

DOUBLE PRECISION DSESD 
RAN1-G6UBFS(DSE£D) 

RAN2=GGUBFS(DSEED) 

GAUSS = (SGP.T (-2»AL0G(RAN1 ) )*C0Sl2*3.14159*RAN2 )*DE)+S 

RETURN 

END 
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C 

SUBFiOUTiNE REACN2 
C 

C FRQCRAK FULheT USES THIS SUBHCUT1N5 TO READ IN IhTECRAl, 

C INTEGRAL ENERRV LdSS DISTRIBUTION FUACIIOWS FOli THE 

C DEPLETION LA^tcd, DEFINED IK GTHIN/ AND EACH SEGMENT 

C OF THE SUBSTRATE TRACK/ DEFINED IN GTHICK. THESE FILES 

C ARE CREATED IN VAV. 

C 

COH»'aN/SIMPLE/(iTHIN(lUC)/GTHICK<lOO)/NOMl/NU»(.2 
COMMON/ YVALUE/ai( IOC)/ TL2 (IDO I 
TYPE 10 

ID FOP.MATdK/lX'SUBROOTINE RSADN2') 

□PtN(UNIT=21/FILS=:'GTHIN.DAT') 

RSAD(21/») NUMi 

READ(21/*) (YLl(I),GTHIN(I),DUM/IsNOMl/l/-l) 
CLaSE(UNIT=2i) “ 

TYPE 20 / NUr'l 

20 FORMAT (IX/IX'GTHIN READ IN ... NPTS * 'lO) 

0PEN(UNIT=:21,FILS='GTHICK.DAT') 

READ(21/*) NUM2 

READ(21/*) < YL2(I)/GTHICK(l),DUM/I=NUM2/l/“l) 

CLOSE (UNIT=21) 

TYPE 30/ NOM2 

30 FQRMATdX/lX'GTHICK READ IN ... NPTS = '13) 

ftSTURN 
END 
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SUBIiOUTINE IdNiZ(9ETA/THICK^ZA/ZX#0EN/DELBAR/SI6,EHAJI/G,BL> 
SUBPOUTINE lONIZ 

j 

PUB POSE 

CALCULATE ICNIZAIIQK LOSS FOR PRCIOKS USIKG BETHP.-BLOCH , 
FORi^ULA WITHOUT DENSITY EFFECT OR SHELL CORfiECTIOh 

DESCRIPTIQB OP PARAHETbRS 
BETA 
THICK 
ZA 
ZI 
DEW 
DSLBAR 
BIG 
EHAX 

G 

BL 

ERM=;511 

THKsI.0E«4*THICK*OEN 
BET A2=BSTA**2 
EKAX=2*ERM*'BETA2/(1-BFTA2) 

GsO .3C‘<'S8»»ERF*ZA’»THK/ (E3TA2^E PAX ) 

FAC=G*EHAX 
ARG=(EBAX/2I)**2 
DSLEARsPAC»(ALQG(AR6) • 2*8ETA2) 

SIGaBHAX^SORT(G*(l - BETA2/2^)) 

BLstDSllBAR*C.6748/FAC»»2 I SILICON 

RETURN 
END 


- VELOCITY OF PROTON /SPEED OF LIGHT 

- thickness in MICBOhS OF ABSORBER 

- ATOMIC NUPbER/ATQPIC MASS OF ABSORBER 

- mean ionization potential 

« DENSITY OF ABSORBER 

- mean ENERGY LOSS (DE/DX) KEV/MICROA 

- STANDARD DEVIATION OF MF.AH ENERGY LOSS 

- HAaIMUH energy TRANSFER TO ELECTRON IN 
SINGLE COLLISION 

- kappa factor of VAVILOV AND LANDAU 
» ELUNCK-LEISEGANG FAC10R 


Oo on 


FUNCTION BXACEHERGY) 
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INPUT IS KINSTIC ENERGY OF PROTON 
OUTPUT IS BETA»V/C. 

PRMa938,213 

ROOTS (PRH/('PRM ♦ SN£RGY))**2 
BTA*SQRTa-ROOT) 

RETURN 

END 


I 
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C 

C PROGRAM VAV.FOR 

C . 

C (iBNERATES DIFFERENTIAL AND INTEGRAL ENERGY LOSS DI STRIBUTir/KS 

C FOR bOTH VAVILOV AND BLUNCK-LBISEGANG THEORIES. PLOTS ON 

C TEKTRONIA AND >»RITBS DISK FILS FOR USE IN PRQGEAri POLHGT 

C 

DIMENSION xaO:),XV(10C>/AKSV(lC:)/YaOC),YSUNaC‘C) 

DIMENSION YH(1C“5),YRSUM(10?) 

DIMENSION LABU59)/LAB2(59),LAB3(59)/LAB4(6),LAH5(8) 

DIMENSION LAB0(59),LAB6(5) 

COMMGN/STTJFF/NCOMP/CCMt/f)Y/NMAX,»-AXB,MAXT/HAXY 
CQHMON/H£TKAP/BSU/D,PVyPf"IN 

COMMON/PLQX/ N/ XMN, XMX, YMN/ VMX/ XI/ X2, X3/ X4/N P/ N C/ 
ILT/ISTP/ISYM/ITX/ITY/IG 
DOUBLE PRECISION FINANE 
C 

C PLOTTING LABELS 

C 

DATA' LABO/'P*/ *r'/ 'o'/ 'b'/ 'a' / 'b'/ 'i'/ '1'/ 'i'/ 't'/ 'y '/ 

148 *' '/ 

DATA L A B3/ 'E '/ 'h '/ *6 '/ 'r */ *9 ' / 'y ' / ' '/ 'L'/'o'/ 's'/'s '/ 

1' *’/'('/'k'/'e'/'v'/')'/42*' '/ 

data L a B1 / * L */ * a */ * 11 ( 1 ' / ' o */ ' q' / *8 * / ' '/'(*/ *L'/ 'a'/ 'n'/'d'/ 
I'a'/'U'/ ')'/44*' '/ 

DATA L AB^/ 'L '/ ' a '/ * n '/ 'o'/ ' d' / ' a ' / ' */*I */*V*/'a'/'v'/'i'/ 
l'l'/'o'’/'v'/')'/43*' '/ 

DATA l.AB4/'K'/'a'/'’p'/'p'/'a'/'='/ 

DATA LABS/'B'/ 'e'/ 't '/ 'a'/ **',*»'/ *2'/ 

DATA LABb/'B'/'*'/'»'/'2','='/ 

C 

C INITIALIZE PLOTTING PARAMETERS 

C 

X1=1.;X2=1.;X3=7.75;X4s5.75;MP='YES';NG='YSS*;LT=C;ISTP=1; 
1I3YM=0;IG=,2.; IXX=1;ITY=1 

c 

C INITIALIZE PARAMENTERS FOR 3-PQlNT SIMPSON RULE INTEGRALS 

C 

NCOMPslC 
C0MP=1.D0S-C5 
NMAX=1()01 
NAXy>C 
MAXBsO 
C 

C PROVIDE INPUT OPTION 

C 

TYPE 5 

5 FORMATCIX/IX'INPUT OPTION! ENERGY-THICKNESS (1) OR KAPPA- 

1BSTA»*2 (2)!'$) 

ACCEPT »/ OPTION 
IF(OPTION.EQ.l) GO TO 15 
C 

C INPUT KAPPA AND bETA**2 

C 

TYPE 25 _ 

25 FORMATdxnX'ENTER KAPPA AND BETA**2!'$) 

ACCEPT */ D, BSQ 
GO TO 3U 
C 

C INPUT PROTON ENERGY AND SILICON ’iAFLfi THICKNESS 
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IS 

lb 


20 

C 

c 

c- 


30 


C 

c 

c 

150 


C 

C 

C 

35 

C 

C 

C 

C 

C 

c 


100 

c 

c 

c 

105 
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TYPE 1C 

fORVATdX/lX'ENTER PROTON KINETIC ENERGY (MEV):'S) 
AO’CSPT *, ENERGY - 
TYPE 20 

FORNATUX/IX'EuTER silicon thickness (AICR3NS):'$) 
ACCEPT »/ THICK 

DEFINE SOME HOPE VARIAHLSS NEEDED IN INTEGRATICN 


Zsl4;As28,0 85;2A=Z/A;D£NSE=2.328AZI=.nC I SILICON 

£>ETA=BTA(ENERGY) 

CALL I0HlZ(BETA/THICK/2A,ZI,DENSE^0SLBAR,DUH/hy0/BL) 

BSQ=BETA»*2 

ERM=511 

»a2*ERM*BSa/(l-BSQ) 

FAC1=D*W 

FAC2=1+BSQ-0.5772164AL0G(DJ 

DY=0.C2/D 

SUM=0«0 

ESTABLISH LiKllS OF CONFUTATION 
TYPE 150 

FORMATCIX/IX’ENTER PMIh/PMAX: *S) 

ACCEPT PMIN/PNAX 
Ns2’»(PMAX - PMIN) 1 

MHAT TO COf<PUTE 

TYPE 35 

FOHMATUX/IX'COMPUTE DIFF (1 ), INT (2 J CR 80TH(3) ?S'$) 
ACCEPT */ LP 

accumulate POIliTS FOR GRAPH 

IF(LP.EQ.2) GO ID 1C5 

DIFFERENTIAL SFSCTRUM 

DO IOC 1=1 /N 
PsPMlN + (I-1)*C.5 
Y(I)=VAVDIF(P) 

YBa)=BLDIF(P,BL) 

X(I)=P 

XV(I)=PV 

XKEV(I)=DELBAR + FACi*(P+FAC2 ) 

CONTINUE 

IF(LP.SQ.l) GO TO 25‘-' 

INTEGRAL SPECTRUM 

MAXY=C 

MAXB=D 

MAXTss* 

DO 110 I=UN 
P=PFIN + (I-l)*0.5 
YSUf(I)sl. - VAVIKT(P) 

YfiSUM(I)=l. - flLINT(PMIN/P/BL) 

X(I)=P 
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XVOsPV 

XKEV(I)=D5LBAP + FACl*(P*^FAG2) 
no CONTINUE 

25C TYPE 30 :> 

30C FORMATUa/IX'X values: landau (1), VAV!'0V(2) CR K2V(3):'$) 

ACCEPT *, NVAL 
TYPE 35C 

350 FORNATdX/lX'EATER X«IN/XMAX/ YMIU/YHAX') 

ACCEPT XMN/XMX/YNN^YMX 
IF(NVAL.^'^.l) CALL PLTGEN(X,Y) 

NP='N0';NG='Na' 

IF(LP.Ea.2.0Pt.LP.Eg.3*AND.NVAL.SC.l) CALL PLTCENCX, YSUl' ) 

C IF(NVAL.EQ,1) CALL PLTGEK(X/YB) 

IF(LP,EU.2.0R.LP.EU.3.AND.NVAL.E(;.1) CALL PLTGENC X/ YHSUN) 
IF(NVAL.SQ-l) CALL LABL(3./0. 25, 'H'/LAB1,0. ,1, 6,1 .5/10/0 ) 

NP= 'YES ';N&= 'YES' 
iF(NVAL.E«.2) CALL PLTGEN(XV,Y) 

NPs'N0';ng=*’K0' 

1F(LP.£Q.2.0R.LP.£Q«3.A.Y0.NVAL.E0«2) CALL PLTGEIKXV/ YSUN) 
IF(NVAL.Ea.2 J CALL PLTGSH (XV/ Yb ) 

IF(NVAL.FQ.2) CALL LABL(3./0. 25/ 'H '/LAB2/0. /I. 5,1 .5/ lC/0 ) 
NP='YES';NGs'YES' 

IF(KVAL.EQ.3) CALL PLTGEN(XKEV/Y ) 

NP^'NO'’;NG='NO' 

IF(LP.EQ.2.0R.LP.E0.3.AND.NVAL.EG.3) CALL PLTGEN( XKE V, YSUM) 
IF(KVAL.BU.3) call PLTGSN(XKI,V,YB) 

IF<LP.EQ.2«0R.LP.Ba.3.AND.NVAL.EU.3) CALL PLTG£N(XKEV/ Y6SUH) 
IF(NVAL.EQ.3) CALL LABL (3./0. 25/ 'H'/LAfi3/0. /I . 5/1. 5/ 10,0 ) 
CALL LABL(0.25,2.75,'V',LAB0,0., 1.5, 1,5, 10,0) 

CALL LAPL(l.,5.8,'H',LAe4,D,1.5,1.5,10,0) 

C CALL LABL(«l.5,5,8,'H',LAb5,BSO,1.5,l,5,10,? ) 

CALL LABL(4.5,5.8, 'H',LAB6,3L,1.5,1.5,1?,0> 

NP='YES*;NG='YES' 

CALL TINPUT(KK) 

TYPE 400 

400 FORNATCIX/IX'ENTER 1 TO GRAPH AGAIN: '$) 

ACCEPT *, GRAG 
IF(GRAG.SO.l) 60 TO 250 
TYPE 50C 

50C FORNATCIX/IX'NAME OF OUTPUT FILEt'^) 

ACCEPT 41, FINAME 
41 FORMAT (AlO) 

OPEN(aNIT=21,FILfc=FlNAM£) 

WRITE(21,*> N 

WRITEC 21,*) (X (I ), YBSUM (I ), YB (I ), 1=1, N ) 
CL0SE(0NIT=21,HLB=FINAMP) 

STOP 

END 
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FOMCIION VAVINT(P) 

C 

C INPUT IS VAVILOV LAMUDA 

C aUTPUT IS 3XPRESSI0K (III.l) IN APPENDIX III 

DIMENSION FIHIOIC ),F2F(1310),FACl(1010)/FAC2(10lCf>/ 
1 FAC3( 1010)/ GO J/ SI U'?1C )/Y(1010 ) 

COKPON/STUFF/NCCNP/CCMP/OY/NMAX/*tAXB,MAXT/MAXY 

COMMON/BETKAP/eSU,D/PV/PMIN 

Y(1)=0.0 

DL=ALOG(D) 

PV=D*(P + DL) 

PVMIN=D*(PMIN + DL) 

GOHG-E.XPC-.0»a.C ♦ 0.5772157*BSQ)> 
CsDY/(3.0*3,141593*G0RG) 

SsQ.C 

NCQONT^O 

G(l)=GORG*(PV - PVMIN) 

.DO 26C L=2/NMAX 
IF (HAXY-L) 190/200,200 
190 Y(L)=Y(L-1) + DV 

CALL TRIGEX(Y(L)/3I(L)/CI) 

FAC1(L)=AL0G(Y(L)) - Cl 
FAC2(L)=C0S(Y(L)) f Y(L)*SI(L) 

FAC3(L)=Y(L)-*FAC1(L) + SIN<Y(D) 

200 IF (MAXB-L) 21C/220/220 

210 FlF(L)sSXP(D'UESa*FACl(L) - FAC2(L))) 

F2f (L)-D*(FAC3(L) + BSQ»SI(L)) 

220 ND*L - 2^(L/2) 42 

SIN1=:SIN(Y(L)*PV 4 F2F(L>) 

SIN2=SIN(Y(L)*PV«IN 4 f2F(D) 

G(ND)=:F1F(L)*(SIN1 - SIN2)/Y(L) 

GO TO (260/260/230 ,ND 
230 ADDsG(l) 4 4.0*G<2) 4 G(3) 

G<1)=G(3) 

SsS 4 ADD 

IF (ABS(ADD/S)>CQMP) 240/250/250 
24C NCOONTsNCOUNT ♦ 1 

IF (NCOUNT-NCOF.P) 260/270/270 
25C NCOUNTsO 

260 CaNlINUE 

MAXT=NHAX 
GO TO 280 
270 MAXT=L 

280 IF (MAXT-MAXY) 300/300/290 

29C .MAXY=MAXT 

300 IF (MAXT-MAXB) 320/320/310 

310 MAXB=MAXT 

32C SUH=C*S 

VAVINT=SUM 

RETURN 

END 
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C 

C INPUT IS VAVILOV LAPBDA 

C OUTPUT IS EXPRESSION DEFINED IN APPENDIX II 

C 

DINENSIOf! Flt-(ini0)/<?2K101C),PACl(lClC),FAC2(lC10)^ 
1 FAC3(lCl'l)/G(3)/SI(101C)/Y(lCilO) 

COMMUN/STUFF/NCOHP/CUMP/DY/NMAX/NAXB/MAXT^MAXY 
CaMMON/BETKAP/BStl/D/PV/PMIN ’ 

Y(1)=C.C 

0L=AL0G'!:D) 

PV=D»(P + DL) 

CORG=EXP(-D*(1.0 + O.S772157»esQ)) 
C=D*DY/<3.C*3.141593*G0RG) 

S=0,0 
NCOUNT=0 
0(1 )-6QR0 
DO 260 L-2/NMAX 
IF (MAXY-L) 190/200/200 
19C Y(L)=Y(L-1) + DY 

CALL TRIGEX(Y(L)/SI(L)/CI) 

FACKD-ALOGCYCD) - Cl 
FAC2(L)=COS( Y(L)) + X(L)»SI(L) . 

FAC3(L)=Y(L)*FAC1(L) + SIN(V(D) 

200 IF (MAXB-U 210/220^220 

210 F1F(L)=5XP(D*(ESQ*PAC1(L) - FAC2(L))) 

F2F(L)=D*(FAC3(L) + flSQ*SI(L)) 

220 ND=L - 2*(L/2) 42 

G(ND)fF1HL)*C0S(Y(L)*PV 4 F2F(D) 

GO to' (260/260/230) /ND 

230 ADDrG(l) 4 4.0«G(2) 4 G(3) 

G(1)=G(3) 

S»S 4 ADD 

IF (ABS(ADD/S)-COMP) 240/250/250 
240 MCQUNT=NCOUNI 4 1 

IF (NCOUNT-NCOHP) 260/27C/270 
250 NCOUNT=0 

260 CONTINUE 

MAXT=NMAX 
GO TO 280 
279 MAXTsL 

28C IF (MAXT-MAXY) 390/3C-0/290 

3QG (4AXV-MAXT 

30C IF (MAXT-MAXfi) 320/320/310 

310 HAXPsMAXT 

320 SUM=C4S 

VAVDIFsSUM 

RETURN 

END 

C 

C 

C 

SUBROUTINE 7RIC-EX(X/SI/CI ). 

X2 = X*X 

IF (X2 -100.0) 10/40/40 
Iw SI=v.O 

CI=C. 57721566 4 G.5*AL0G(X2) 

T-X 

Asl.O 
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20 SI=SI * T 

AsA l.Q 

I=-T*(X’»(A.-1.0)/(A^*2)) 

IF (ABS(T)-(1.0S-»>9)) 100/30/30 
3C CIS Cl ♦ T 

AsA 1.0 . 

T=T*(X*(A-1.0)/(A*»2)) 

IF (AB?(T)-a.CS-0y)) lCO/20/20 
40 AMAXs0.5»(X ^ l.v) 

Psl.C 

QsO.Q 

Ssl.O/X 

Asi.O 

50 Us{| s 

A=A/^*^.C- 

IF (A-AMAX) 60/60/90 

60 S=-S*(A/X) 

P=P+S 

IF (ABS(S)-(l.uE-12)) 90/70/70 
70 A=A 4- 1..J 

IF (A-AMAX) 80/80/90 

80 SsS*(A/X) 

GO TO 50 

9C CXsCOS(X)/X 

SX=SIN(X)/X 

SIsl. 5707963 - P»CX - Q*SX 
ci=p»sx - a^cx 
100 RETURN 

END 
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HfNCTIQJI BLDIP(P/Bi) 

D I M e t-S ION C ( y )> SL'-'BD A ( 9 ) ^ fi AMM A < 9 ) 

L’t A T A C / . • 3 1) a / . •: 9 -1 3 / . i, 8 8 2 / , : 6 ^ 7 / . C 3 5 9 , . . 1 6 4 , * C 6 4 / . C (• 2 1 , 
Df.rf RLMBDA/-! .46/-:.738y .17/ 1.3a,2.S5/5.39/S.4/if .8/3% 
iVATA GA''‘iA/. 737/. 947/1. 23/1.08/2. 4/3.68/6.18/12.3/30.7/ 
SU.V-'..; 
r^. n 1 1 s* 

KAC1=C(I)*GAMMA(2) 
bAC2sGAK*-!A(I I**? ♦ P2 
F?C3=P - RLMBD<t(l) 

SUH^aUM ♦ FACl»EX?(-FAC3*«2/FAC2)/SChT(FAC2) 

c'intivsj: 

HLDlFsS'JiN 
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Figure 1. Cutauay view of diodes. P+ layer la about 0.7 «. Vertical scale la eiaggetated. 










Figure 3. Test configuration used 
rotated In 0.1“ incremei 
or out of the proton be; 




GeooietricsX axiTdngeoient of diodos on t^lio axxay. 
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Figure 9. A Z-intensif led contour of pulse height of coincident pulses in 
one diode versus pulse height of poincldent pulses in a second 
diode. 9(a) is for a pair of diodes separated by an even number 
of diodes while 9(b) is for a pair separated by an odd number of 
diodes. 
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Diffusion Lsngtb 

Figure 12. Charge collection efficiency vs. diffusion length for 400 y track entering at x » y *0. 

Wafer thickness 400 y. (a) 0-100 y (b) 100-200 y (c) 200-300 y (d) 300-400 y. °Soll2 Unl 
— V — 0, Dashed line -- v « 
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Figure 14, Diff 
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Figure 15. Diffusion current with 
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Figure 16. Diffusion current with 
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Figure 17. Diffusion current with L 
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Figure 18. Diffusion durrent with I», = 25 p 
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0.50 1.75 3-00 4-25 5.50 6.75 8.00 9.25 10-50 11.75 

ChorgR Coffected CfC) 

.Figure 24. Computed Integral pulse height spectrum for 21.5 McV protons using Gaussian f(A,s) along 
400 u track. 


ORIGINAL PAGE IS 
OF POOR OUAUTY 






ORIGINAL PAGE IS 
OF POOR QUALfTY 



s^unoQ 


98 




33 3 




Figure 30. Protoa entry points corresponding to charge 3.2 fC collected by Diode at 21-5 MeV. 
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